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ABSTRACT 

We present the results of a multi-wavelength multi-epoch survey of five evolved protoplanetary disks 
in the IC 348 cluster that show significant infrared variability. Using 3-8/xm and 24/im photometry 
along with 5-40/xm spectroscopy from the Spitzer Space Telescope, as well as ground-based 0.8-5^m 
spectroscopy, optical spectroscopy and near-infrared photometry, covering timescales of days to years, 
we examine the variability in the disk, stellar and accretion flux. We find substantial variations (10- 
60%) at all infrared wavelengths on timescales of weeks to months for all of these young stellar objects. 
This behavior is not unique when compared to other cluster members and is consistent with changes 
in the structure of the inner disk, most likely scale height fluctuations on a dynamical timescale. 
Previous observations, along with our near-infrared photometry, indicate that the stellar fluxes are 
relatively constant; stellar variability does not appear to drive the large changes in the infrared fluxes. 
Based on our near-infrared spectroscopy of the Pa/3 and Br7 lines wc find that the accretion rates are 
variable in most of the evolved disks but the overall rates are probably too small to cause the infrared 
variability. We discuss other possible physical causes for the variability, including the infiuence of a 
companion, magnetic fields threading the disk, and X-ray flares. 



1. INTRODUCTION 

The general picture of the formation and evolution 
of low mass stars involves the complex interplay of the 
newly created star and the reservoir of dust and gas 
surrounding it, feeding it and often shielding it from 
out view. Early on, material from a dense envelope 
surrounds the protostars, falling onto the disk from 
where it accretes onto the protostar building up a ma- 
jority of the stellar mass (|Terebev. Shu fc CassenllTOSl 
I Adams. Lada fc Shulll987f ). After hundreds of thousands 
of years the envelope dissipates, the protostar has ac- 
creted most of its mass, and the system is left with a disk 
of gas and dust. These circumstellar disks are usually vis- 
ible for millions of years before much of the dust and gas 
is either accumulated into planetesimal s and planets, or 
remo ved by a photo-evaporative flow (|Hernandez et all 
[200l . Left behind is a debris disk, with a small trace of 
optically thin dust sustained by the collisions of planetes- 
imals (Wyatt 2008). While these different evolutionary 
stages of low mass star formation have been fairly well 
established, the physical processes behind the movement 
from one stage to the next are not well understood. 

Transition disks, and evolved disks in general, repre- 
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sent a stage of the disk's life that is believed to lie be- 
tween fully optically thick disks and debris disks. These 
systems are characterized by a deviation from the typ- 
ical SED of a full flared disk, often due to the pres- 
ence of large, many AU wide regions of optically thin 
material or a severe flattening of an otherwise flared 
disl0. These SED sha pes may be indicati ve of th e pres- 
ence of planets fBr vden et al.| [l999: Nai ita et all [2007). 
a photocvaporative fl ow combined with viscous evolu- 
tion (Alexander et all l2006t iPascucci et al.l l2011h or a 
significant amount of small dust settling to the midplane 
(jCieza et al.ll2010l ). Resolved observations confirm the 
presence o f gaps filled with opti cally thin dust spanning 
many AU ([Andrews et al.|[2011[ ). as well as the presence 
of optically thick material ve ry close to the s tar inte- 
rior to the optically thin gap (jPott et al.ll2010t ) in some 
of these evolved d isks. Gas has been obseryed in many 
of th ese systems (|Pontoppidan et all 120081 : iSalvk et "all 
1201 If ) at levels higher than expected given the drop in 
dust surface density needed to explain the SED. Undcr- 

* Different studies apply different definitions to the term 'tran- 
sition disk' in an attempt to best encapsulate the evolution of cir- 
cumstellar disks. Some authors restrict this analysis to sources 
that contai n a gap i n the disk that is nearly devoid of dust 
llMuzerolle e tTalllSoToD . These sources are characterized by an in- 
frared SED that has little to no excess at wavelengths less than 
^ 10/jm with a strong excess at longer wavelengths. Some authors 
allow for some optically thick material at the very inner edge of the 
gap, close to the dust sublimatio n radius, w hich leads to a strong 
excess in the near-infrared (EsDaillat et al.' 2007). Other studies 
include disks that have merely been flattened, rather than having 
a gap (Cicza et al. 2010). These full disks do not show a sharp dip 
in the infrared SED, but instead show an overall steepening until 
it almost resembles the emission from an infinitesimally thin, per- 
fectly flat disk. All of these disk geometries represent some form 
of dust evolution, but not necessarily the same type of evolution. 
Here we choose to use the term 'evolved disks' to describe any 
system that shows significant evolution in its dust content, while 
'transition disk' only describes those systems with either a hole or 
gap. 
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standing the gas and dust in these situations is crucial for 
a complete picture of planet formation and early stellar 
evolution. 

Recently some transition disks have been found to 
have large amplitude infrared variability, which could 
potentially be relate d to the clearing of the disk. 
iMuzeroUe et all H20Q9f ) study one member of the IC 348 
cluster in which the flux shortward of 8/im increases 
while the flux longward of 8/im decre ases (and vice- 
yersa), both in as little as one week. lEs paillat et al.l 
see this 'seesaw' behavior in a collection of 14 
transition disks in Taurus and Chameleon observed with 
Spitzer 5-40/im IRS spectroscopy. The strength, speed, 
and wavelength dependence of the variability was sur- 
prising given the previously known sources of variability 
in pre-main sequence stars, such as starspots and stel- 
lar occultations. Even large outbursts, such as seen in 
FU O ri or EXor objects, take place over mo nths and 
years HHartmann fc KenvonllTOQa lHerbisj[2008l ). Fluctu- 
ations in the 20/xm flux on weekly timescales are unex- 
pected given that this wavelength traces dust at a few 
AU where the dynamical timescale is closer to years. 
Recent observations suggest that rapid infrared variabil- 
ity may in fa ct be quite common arnong p re-main se- 
quence stars ([Morales- Calderon et al.l [2011D . Geomet- 
ric models of the disk suggest that this variability may 
be due to changes in the scale heig ht of the inne r disk 
(iFlahertv fc Muzerolle 2010; Espaill at et 3111201 ID 

In ' Flahertv et al.l (|201 1) we presented the results of a 
detailed observing campaign focused on the source LRLL 
31 in IC 348. Using infrared observations from Spitzer, 
as well as ground-based near-infrared spectroscopy, we 
were able to confirm the geometric model of inner disk 
perturbations, as well as rule out accretion as the driv- 
ing source of the variability. Here we present multi- 
wavelength multi-epoch observations of five additional 
evolved disks in IC 348 designed to determine if the be- 
havior observed in LRLL 31 is common, as well as put 
further constraints on the source of the vertical pertur- 
bation. These additional sources were chosen to have 
strong infrared variations as well as evidence for signifi- 
cant evolution in their dust structure based on the SED. 
In Section [5] we detail the data obtained for our sample, 
review the properties of LRLL 31, and discuss in detail 
the characteristics of the stellar, dust and gas emission 
for each source. We then use this information in Sec- 
tion [3] to determine if the same process is operating in all 
of these stars. Assuming this process is common among 
normal disks, in Section 2] we discuss what could be its 
physical cause. 

2. VARIABILITY OBSERVATIONS 
2.1. Data 

We obtained multiple epochs (2007-2009) of data from 
4000A out to 40^m on the IC 348 evolved disks LRLL 
2, 21 , 58, 67, 1679 ()Luhman et al.l [20031 [M^ench et al.l 
|2007| ) as part of the same program t hat observed t he IC 
348 transition disk LRLL 31 (Flah ertv et al.l 1201 Ih . We 
also include previous infrared Spitzer observations taken 
as early as 2004 to extend our temporal coverage. Table[T] 
lists all of the observations for each source, including 
the times when these measurements were taken. These 
observations include: 



• 10 epochs of Spitzer 24/zm photometry with typical 
uncertainties less than 2% (Table [2]). Two epochs 
were taken seven months apart in 2004 as part of 
the c2d map of IC 348 (jRebuU et al.ll2007D while 
another epoch wa s taken in 2004 as part of GTO 
observations (Lad a et al.l [200(1 ). We obtained five 
consecutive days of monitoring in Sep 2007 followed 
in Mar 2008 by two epochs separated by one week. 

• 7 epochs of Spitzer IRAC cold-mission 3.6,4.5,5.8 
and S.Ofim photometry with uncertainties less than 
2% (Table [3]). One epoch was ta ken in 2004 as 
part of the c2d coverage ( Jorgens en. et al.l I2006D 
while GTO observatio ns were also performed in 
2004 ([Lada et al.ll2006D . We obtained five consec- 
utive days of photometry in Mar 2009. 

• 4 epochs of Spitzer IRS 5-40/j,m spectra. Two spec- 
tra were taken a week apart in Oct 2007, while two 
more were taken in Mar 2008, also separated by 
one week 

• 38 epochs of Spitzer IRAC warm-mission 3.6 and 
4.5/im photometry spread over 40 days in Fall 2009 
with uncertainties of 1-2% (Table Hj). The cadence 
for this photometry ranged from four observations 
a day to one observation every two days throughout 
the visibility window. 

• One epoch of low-resolution optical spectra taken 
with the CAFOS instrument on the 2.2-m telescope 
of the Calar-Alto Observatory of LRLL 21, 58 and 
1679 in Fan 2009. 

• Multiple epochs of 0.8-5//m spectroscopy with Spex 
on IRTF in FaU 2008 (LRLL 21, 1679) and Fall 
2009 (LRLL 2, 21, 58, 1679). The exact number 
of epochs and the wavelength coverage of the spec- 
tra varies from source to source, with each object 
getting at least two 0.8-2.5/im spectra. 

• One epoch of high-resolution optical spectra of 
LRLL 21, 58, 67 and 1679 taken with Hectochelle 
on the MMT in Feb 2008. 

Detai ls on the dat a reduction can be found in 
Flahe rtv et al.l ([201 If) . We also observed four non- 
accreting diskless weak line T Tauri stars (WTTS) in 
Taurus (KO star LkCa 19, A2 star HD 57928, M3 star 
LkCa21, Ml star JH 108) with Spex to use as standards 
for the 0.8-2.5/zm spectroscopy. We have found that 
pre-main sequence stars provide much better templates 
than the dw arfs and giants in the IRTF spectral library 
([Ravner et al. 2009) when comparing the spectral shape 
and strength of the absorption lines. This is especially 
true for the spectral shapes of the late-type stars and the 
absorption lines in the early type stars. For the IRS spec- 
tra, unlike with LRLL 31, we did not obtain both Short- 
Low (SL) and Long-Low (LL) spectra at every epoch, 
due to the significant background emission that made it 
difficult to extract the LL spectra. For LRLL 2 and 58 
we obtained four epochs of SL spectra, but only included 
LL spectra in the first two epochs. These two sources, es- 
pecially LRLL 58, are the most susceptible to systematic 
uncertainties in the long wavelength end due to the large 
background. The background emission also leads to a 
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systematic offset between the SL and LL spectra that is 
seen in LRLL 2, 58. When comparing spectra taken on 
different epochs we have scaled the LL spectra to match 
the SL spectra. For LRLL 67 the offset maybe due to cal- 
ibration uncertainties causing the short-wavelength end 
of the LL2 spectrum (~ 14 — 16/im) to flatten out. 

We chose to focus on systems with SEDs indicative of a 
significant evolution in the dust distribution from that of 
a full, optically thick, flared disk. Recent radiative trans- 
fer modeling has shown that two of the sources we study 
(LRLL 21 and 31) have large gaps that extend from the 
inner disk out to tens of AU, another source (LRLL 67) 
has no optically thick material within tens of AU of the 
star, while anothe r source (LRLL 2) ha s a flattened full 
disk with no gap (lEspaillat et al.|[2012f) . Given that our 
sources represent a wide range of geometries, we choose 
to use the term 'evolved disks' when describing the en- 
tire sample. Evolved disks can be further subdivided into 
those systems with SEDs characteristic of a flattened disk 
and those systems with SEDs characteristic of an opti- 
cally thin gap, which we refer to as transition disks. The 
transition disks can then be further divided into those 
that only have optically thin material within a few AU, 
such as LRLL 67, and those that have some optically 
thick material at the dust sublimation radius, which we 
refer to as pre-transition disks. As we will discuss in sec- 
tion |3l the most important piece of the disk when it 
comes to the variability is the optically thick dust near 
the sublimation radius, which we refer to as the 'inner 
disk'. 

2.2. LRLL 31 

Here we summarize previously published results o n 
LRLL 31 (jMuzeroUe et all [2009t iFlahertv et all l20TTh . 
The data set is similar to, but more extensive than, the 
other evolved disks. One of the main questions we are 
trying to answer is whether or not the variability ob- 
served in LRLL 31 is common or is unique to this one 
source. We start here with the detailed observations and 
in section [3] we try to determine if we can synthesize this 
information into one simple model. 

2.2.1. Stellar Properties 

LRLL 31 is a G6 star with a lumin osity of A.2>L(r,, 
R^.IRq, M=1.5M0 and Ay = 8.8 (jFlahertv et alJ 
1201 It ). We estimate the extinction using the shape of 
the 0.8-2.5/im spectra, ignoring the K-band whose shape 
may be influenced by dust emission. Contemporaneous 
with the near-infrared spectra, we were able to obtain 
J-band photometry, which is dominated by stellar emis- 
sion, allowing us to calculate the stellar luminosity. Our 
multiple epochs of J-band photometry do not show signif- 
icant fluctuations in the luminosity, although there may 
be long term trends in the stellar flux. The near-infrared 
spectra show that the extinction is constant on timescales 
of weeks and years. There is some evidence for variations 
in the radial velocity of the star, suggesting the presence 
of a companion, although with only three epochs this 
conclusion is not deflnitive. 

2.2.2. Infrared Variability 

Our infrared observations cover A = 2 — 40/.tm on 
timescales of days to years allowing us to look for any 



wavelength dependence in the fluctuations, as well as to 
constrain the variability timescale. In the IRS spectra, 
LRLL 31 displays large fluctuations 60%) in as lit- 
tle as one week with a clear wavelength dependence in 
which the SED appears to pivot around 8/im. The 24/im 
photometry displays large fluctuations (30%) consistent 
with those seen in the IRS spectra. The flux of the sili- 
cate feature also changes with time, and the changes are 
correlated with the long- wavelength fluctuations. The 
3.6 and 4.5/im monitoring shows large (0.3 mag) fluc- 
tuations predominately on weekly timescales. Multiple 
epochs measuring the shape of the 2-5/im excess find that 
this excess looks like a single temperature blackbody at 
~1700K whose strength rapidly varies with time. The 
[3.6]- [4.5] color of the system appears to get redder as it 
gets brighter. There was no evidence for periodicity in 
the light curve. 

To describe the observations, we developed a model 
for disk structural perturb ations to the inner disk 
(jFlahertv fc MuzeroUd IMot) . LRLL 31 is a pre- 
transition disk and the basic geometric picture is of an 
optically thick ring of dust at the dust sublimation ra- 
dius (the inner disk) that emits predominately in the 
near-infrared, a large gap that may or may not be filled 
with optically thin dust, and an optically thick, fiared 
outer disk that dominates the fl ux longward of 10/im 
(Fig [1^). lEspaillat et al.l ()2012D recently constructed 
a detailed radiative transfer model that conflrmed this 
structure. They flnd that the SED can be fit when in- 
cluding a wall of optically thick dust at the sublimation 
radius (0.3 AU), followed by optically thin dust extend- 
ing to 14 AU beyond which is an optically thick disk. 
In our conceptual model of the variability, when there 
is no vertical perturbation to the inner disk (Fig [T}d) its 
flux is at a minimum because its emitting area is the 
smallest, while the long- wavelength flux is at its strongest 
because it is directly heated by the star. When a, possi- 
bly non-axisymmetric, perturbation increases the height 
of the inner disk (Fig [1];) the short-wavelength emission 
from this region increases, because the emitting area of 
the inner disk is larger, while the long-wavelength emis- 
sion decreases, because the regions of the outer disk re- 
sponsible for this emission are being shadowed by the 
inner disk. This model is consistent with the wavelength 
dependence of the variability in the IRS spectra, the 
daily to weekly timescales seen during the warm-mission 
monitoring campaign, and the change in the 2-5/im ex- 
cess, which is domina ted by emission from the inner disk 
(jFlahertv et al.|[20Tll ). The presence of an optically thin 
gap and an optically thick, flared outer disk might not 
be strictly necessary to explain the variability, but we 
include it in the basic geometric picture because it is 
consistent with the general shape of the infrared SED. 

2.2.3. Gas Properties 

Our optical and near-infrared spectroscopy of the hy- 
drogen emission lines (Ha, Pa/3 and Br7), whose flux is 
directly related to the accretion rate, as well as measure- 
ments from the literature, allow us to characterize the 
level of variability in the accretion rate and look for any 
connection with the infrared variability. Based on multi- 
ple epochs of the Pa/3 and Br7 emission lines, we observe 
a factor of five change in the accretion rate, from 0.3 
to 1.6 X lO~®M02/r~^. Accretion rates are derived from 
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Fig. 1. — Schematic diagrams (not to scale) describing the basic 
geometric model for the variability, (a) We consider an equilibrium 
situation that consists of an optically thick, flared, outer disk, an 
extended gap possibly filled with some optically thin dust, and 
optically thick dust near the dust destruction radius. This last bit 
of dust close to the star likely has a small scale height, and may 
not be axisymmetric, in order to match the weak average short 
wavelength excess seen around most of our evolved disks, (b) The 
variability, in its simplest form, involves an oscillation between two 
states. In one state, the inner disk is very flat and the outer disk 
is directly illuminated, (c) When a vertical perturbation exists in 
the inner disk, then its flux increases because its emitting area is 
larger, while at the same time it casts a shadow on the outer disk, 
reducing its temperature and flux. This perturbation may or may 
not be azimuthally symmetric. 

Pa/3 and Br7 EW using the JHK photometry combined 
with the veiling as a n estimate of the continuum level 
(|Flahertv et al.ll201lD and then converti ng from line flux 
to accretion rate (jMuzeroUe et al.lll99"8l ). Accretion rate 
variability is also seen in multiple epochs of Ha spec- 
troscopy, which show large fluctuations in the EW. The 
accretion luminosity is small (Lqcc/L*< 0.1) compared 
to the stellar flux. While there is evidence for a corre- 
lation between the accretion rate, as measured with the 
near-infrared spectra, and the infrared flux, as measured 
by our warm-mission monitoring, in our 2009 data we 
found that the change in accretion luminosity was too 
small to explain the structural perturbation. 



2.3. LRLL 2 
2.3.1. Stellar Properties 

LRLL 2 is a pre-main sequence A2 star with a bolomet- 
ric lu minosity of 137^0, a radius of SR© (iLuhman et al 
l2003f) and a mass of 3.5Mq based on the iSiess et al 
(jlOOO) 3 Myr isochrones. In our two epochs of near- 
infrared spectra we measure an average extinction of 
Ay=2.9 (Table El) , which is c onsistent with previous es- 
timates ([Luhman et al.|[2003D given that we use a differ- 
ent extinction law (Ri/=5.5 vs. 3.1, see iFlahertv et al.l 



((20ll for details). Previous studies have found no 
evidence for variability in the opti cal, suggesting that 
the photospheric fl ux is steady (jCohen et al.l |200"^ 
ILittlefair et all[2nn5h . 

2.3.2. Infrared Variability 

The IRS spectra of LRLL 2 show large 20%) vari- 
ations in as little as one week and the change in the 
SED shows a wavelength dependence where the SED ap- 
pears to pivot at A = 6/im (Fig [2] and [3] summarize all 
of the Spitzer data). There is a systematic offset be- 
tween the LL and SL spectra (^ 10%) and we reduce 
the LL spectra by this amount when comparing spectra 
from different epochs. To determine if the silicate flux 
changes in addition to the continuum we subtract a sec- 
ond order polynomial fit to the continuum (defined as 
5 < A < 7/im and 13 < A < 14/im) from the silicate 
feature to produce the result in Figure [H The flux of the 
silicate emission feature changes between these spectra, 
and is proportional to the long-wavelength continuum, 
similar to LRLL 31. The 24/im photometry shows flux 
variations up to 25%, similar to that seen in the IRS 
spectra. Our five consecutive days of 24/im monitoring 
show little change, suggesting that the largest variations 
are on longer timescales. This is consistent with the five 
consecutive days of 3-8/xm photometry, which show al- 
most no variation. Our 3.6,4.5/im monitoring offers a 
more complete picture of the timescale and shows slow 
variations (~ 0.1 mag) on weekly timescales with no de- 
tectable periodicity. 

These long-wavelength data trace material many AU 
from the star while our ground-based spectroscopy is sen- 
sitive to much warmer dust close to the star, assuming 
thermal equilibrium temperatures. Based on the shape 
of the 0.8-2.5/im spectra, there appears to be excess emis- 
sion in K band (Fig [4|) . We estimate the veiling to be 
^ 0.2 in the K-band based on the shape of the spec- 
tra0. This is weaker than is typically seen around stars 
with dust in a puffed inner rim of a disk. We do not 
have enough spectral information to determine if the ex- 
cess is due to hot dust at the dust destruction radius, 
or optically thin material further out, but we can say 
that the excess, including the 3.6 and 4.5/Ltm photome- 
try, is consistent with optically thick dust at T>1000K. 
HspaiUat et al. (2012) find that the infrared SED can be 
fit with optically thick material that extending inward to 
the dust destruction radius (1.68 AU). There model does 
not require an optically thin gap and instead the disk ex- 
tends continuously out from the dust destruction radius 
and the diminished infrared flux relative to a normal disk 
is due to signiflcant dust settling. 

2.3.3. Gas Properties 

During two epochs of observations separated by 3 
weeks in 2009 we see little change in the strength of 
the Pa/3 and Br7 lines (Fig [5]) . The line strengths are 
weaker than simple photospheric absorption, suggesting 
additional emission within the system, most likely due 
to hot gas in the accretion flow. In Figure [5] we dis- 
play the lines after the photosphere has been subtracted, 

^ Veiling cannot be determined from line EWs because of the 
absence of non-hydrogen photospheric absorption lines in the A2 
stellar spectrum 
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Wovelength (^im) Wovelength (jirvi) Wavelength {/im) 

Fig. 2. — Spitzer spectra of LRLL 2. Left: The spectral energy distribution showing the IRS data (colored lines), the IRAC cold mission 
data (dots), JHK photometry (triangles), optical photometry ( squares) and the range of the IRAC warm-mission photometry (error bars 
at 3.6,4.5/xm). A reddened stellar photosphere from .Kenvon fc HartmannI l|1995 ) is shown as the dashed line for comparison. Middle: 
Difference spectra between different epochs of IRS spectra. Black line is the change from epoch 1 to epoch 2 (separated by one week), red 
line is the change from epoch 2 to epoch 3 (separated by 5 months) and the blue line is the change from epochs 3 to epoch 4 (separated 
by one week). Right: Residual flux in the silicate feature after the continuum has been subtracted. 
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Fig. 3. — Spitzer photometry for LRLL 2. Top left Cold mission light curves from our five consecutive days of monitoring along with 
earlier observations (left of dotted line). Different wavelengths are marked by different colors and lines, with the flux normalized to the 
median in each band. Top right: Warm-mission light curve. The error bars show the uncertainty in the 3.6 and 4.5/im warm-mission 
photometry (left and right error bar respectively). Vertical bars show the dates when near-infrared spectra were obtained. Points to the left 
of the dashed line arc from our five consecutive days of monitoring, as well as previous observations during the Spitzer cold-mission. Bottom 
left: MIPS 24/im light curve. Inset zooms in on the five consecutive days of photometry. Bottom right Color-magnitude diagram with 
warm-mission photometry marked with squares, and cold-mission photometry marked with triangles. The error bar shows the uncertainty 
in the photometry. 



showing the emission more clearly, while Table [5] lists the 
measured equivalent widths (EW). The difference in EW 
between the two epochs is most likely due a difference 
in our determination of the continuum rather than real 
fluctuations in the accretion rate. A slight mismatch in 
spectral type can have a significant effect on our subtrac- 
tion of the line, especially in the line wings, but does not 
strongly affect the conclusion that there is some emis- 
sion in the line due to an accretion flow. Accretion rates 
derived from these infrared lines are listed in Table [T) Re- 
cently Donehcw & Brittain (201l[) have shown that the 
relationship between Br7 line flux and accretion rate de- 



rived from low-mass stars, which is used here, is also valid 
for higher mass Herbig Ae/Be stars. The average accre- 
tion rate is 3xlO~^M0yr~^ and the accretion luminosity 
is small relative to the stellar luminosity (Lcicc/L*=0.06). 
The uncertainty in the accretion rate is a factor of 2 due 
to possible spectral type mismatch, uncertainties in the 
location of the continuum and the uncertainty in the con- 
version from line flux to accretion rate. 

2.4. LRLL 21 
2.4.1. Stellar Properties 
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Fig. 4. — Spex spectra for LRLL 2. Spectra have been scaled to 2MASS photometry. IRAC warm-mission photometry has been included 
where available. A stellar photosphere (HD 57928 + a Kurucz model) shown in red has been reddened by Ay=2.9 and scaled to the spectra 
at J band for comparison. Blue points are the photometry, and the error bars in the X-axis show the size of the passbands. 
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Fig. 5. — Pa/? and Br7 lines for LRLL 2. On the left the two epochs (dashed and dotted lines) are directly compared, along with a 
WTTS standard (solid line). There is little evidence for variability between these two epochs. On the right the spectra (solid lines) are 
shown after the photospheric absorption has been subtracted (dotted lines). Once the strong intrinsic absorption has been removed, there 
is evidence for emission indicating ongoing accretion. 



LRLL 21 i s a KO star originally classified by 
iLuhman et al.l (120031 ). whose stellar flux varies substan- 
tially. We were able to obtain seven epochs of near- 
infrared photometry concurrent with our ground-based 
infrared spectra, which includes three epochs in 2008 
separated by one week and 4 epochs in 2009 spread out 
over one month. Details on the reduction of these data, 
along with how the luminosit y is derived from the J band 
photometry, can be found in Flaher ty et al.l ()2011[ ) . We 
also utilize 2MASS photometry taken in 1998. We find 
large year-to-year variations in the J band flux of 0.4 
mag (Table [8|) but no daily or weekly fluctuations larger 
than 0.1 mag. In 2009, when we have contemporaneous 
near-infrared and 3.6,4.5/im photometry, we find no sig- 
nificant change in the J band flux from Oct 8 to Nov 8 
despite a large increase in the mid-infrare d flux over this 
same period. Optic al photometry from |Luhman et ahl 
(|2003h (1=13.21) and lCieza fc Balibed ([2006!) (R=14.69) 
is consist ent with the low stellar flux state seen in 
2008. Cie za k Balibeii ()2006[) find smaU short-term op- 
tical fluctuations likely due to cool spots rotating across 
the surface of the star and deduce a rotation period of 
2.5 days. This is consistent with the lack of large daily 



and weekly changes in the stellar flux in our near-infrared 
photometry. Our estimates of the luminosity from the J 
band fluxes range over 2.52-3.66Lq (Table[|]), which cor- 
responds to M-1.8Mq based on the iSiess et all (j200"0l ) 
isochrones. We measure a constant extinction (Ay=4.2, 
Table [5]) during each epoch of near infrared photometry, 
suggesting that the obscuration is not variable, unless 
the obscuring source has large enough grains that it does 
not redden the central star, although the presence of very 
large grains in the disk is inconsistent with the strong sil- 
icate emission feature. 

There is also evidence that this system is a binary. 
Based on our high-resolution optical spectra, we find 
two peaks in the cross-correlation function, one at -15.7 
km/sec and one at 66.8 km/sec ( Table [TU)l. Dahm (200^ 
measure a radial velocity of -1-17.01 km/sec, but do not 
mention any sign of binarity. It is possible that iDahml 
(|2008f l observed the system during a phase when the two 
components were aligned along the line of sight and there 
was no significant difference in their radial velocities. 



2.4.2. Infrared Variability 
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Fig. 7.— Same as Figure |3] but for LRLL 21. 

In the IRS spectra, LRLL 21 displays large variations 
(~ 40%) at A < 15/ini in observations separated by one 
week with no fluctuations (< 8%) at longer wavelengths 
(Fig [5]). However, the 24/L(m photometry does show fluc- 
tuations up to 20% (Fig [71) , and this flux is systemat- 
ically lower than the IRS spectra, suggesting that the 
long wavelength flux does vary. The flux of the silicate 
feature does not change with time, despite the change in 
the underlying continuum. Our five consecutive days of 
3-8/L(m photometry, as well as our 3.6,4.5/im monitoring 
show rapid variations even from one day to the next, al- 
though the predominate timescale appears to be weeks 
and months with no sign of periodicity. Of all the stars 
in our sample, this one shows the largest fluctuations. 

For LRLL 21 we have ground-based spectra (Fig [5]) al- 
lowing us to measure the shape of the excess emission 
from 0.8-5/im, which in most T Tauri stars is dominated 
by emission from the op tically thick inner wall at the 
dust destruction radius (Muze rolle et al.l I2003a[ ) . Here 
we briefly summarize our proced ure for deriving the ex - 
cess spectrum with more details in lFlahertv et all (|2011l ). 
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First we measure the veiling (r = Fexcess / Fphotosphere) 
by comparing the LRLL 21 spectra to that of a WTTS 
of the same spectral type (LkCa 19) in fifteen small bins 
from 0.8 to 2.5/im. We fit a line to these veiling measure- 
ments as a function of wavelength to derive the K-band 
(2.15/im) veiling. We use this value to normalize the 
LRLL 21 spectra to the photosphere level, whose shape 
is estimated using our KO WTTS standard along with 
a Kurucz model extension (Te//=5250, log g =2.5) be- 
yond 2.5/um. The LRLL 21 spectra and the photosphere 
are then subtracted to produce a spectrum of the in- 
frared excess. The 0.8-5/im excess emission, shown in 
Figures 191101 is consistent with an optically thick disk 
with a temperature of 1900 K at every epoch (Table [TT]), 
while the strength of the emission rapidly varies with 
K-band veiling measurements ranging from 0.19 to 0.56 
(Table fTT]) . The typical uncertainties in the derived 
temperature and veil ing are 200K and 0.1 respectively. 
lEspaillat et al.l ()2012D are able to fit the infrared SED 
with a disk model that includes an optically thick inner 
disk ring at the sublimation radius (0.13 AU), consistent 
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Fig. 8. — Spex spectra (black line) and photometry (blue points) 
for LRLL 21. Spectra have been scaled to match photometry. A 
stellar photosphere (LkCa 19 + Kurucz model) shown in red has 
been reddened by the Ay measured on each night to LRLL 21 and 
scaled to the short wavelength flux of the spectra. Spectra have 
been scaled by a constant factor for clarity. 

with our near-infrared spectra, followed by a gap of op- 
tically thin dust that extends to 9 AU. In 2009 we see a 
large increase in veiling from Oct 8 to Oct 31, which is 
consistent with the large increase in infrared flux seen in 
the 3.6 and 4.5/im photometry during the same time pe- 
riod. The dust temperature we derive is slightly higher 
than expected for the sublimation of silicate dust, which 
occurs closer to 1500 K, similar to the high te mperature 
we observed in LRLL 31 (jFlahertv et al.ll201lD . Not tak- 
ing account of a possible massive companion (discussed 
above) could lead us to overestimate the dust tempera- 
ture by '^100 K and underestimate the K-band veiling by 
'^O.l. It is also possible that the inner disk is not entirely 
optically thick. If the inner disk is azimuthally asymmet- 
ric, with part of the disk being optically thick while part 
is optically thin, then trying to fit this emission as an op- 
tically thick blackbody will lead to an overcstimation of 
the temperature by a few hundred kelvin. Since the tem- 
perature stays constant, these systematic effects are also 
likely constant and should not effect our interpretation 
of the variability. 

The dominant change in the system is in the strength 
of the disk emission, not the temperature of the dust. 
The [3. 6] -[4. 5] color of the system gets redder as the flux 
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Fig. 9. — The infrared excess of LRLL 21. Each plot shows 
the difference between LRLL 21 and the standard, normalized to 
the photospheric flux at 2.15/im. Some have been moved upward 
for clarity and the horizontal dotted line shows the zero point for 
each excess spectra The best fit blackbody is shown with a red 
dashed line and the 3.6,4.5^m photometry is included when avail- 
able. Each spectrum has been smoothed by a median filter 0.01/im 
wide to reduce the noise in the continuum. The parts of the spec- 
tra marked in blue are strongly affected by the telluric correction. 
Note that the scale changes for each plot. The strength of the 
excess shows large variations from night to night, although the 
temperature stays roughly constant. 

increases (Fig[7|). This is consistent with blackbody emis- 
sion from dust, on top of the stellar continuum emission, 
increasing in strength while not changing in tempera- 
ture. Despite being higher than expected, the fact that 
the dust temperature is close to the sublimation tem- 
perature suggests that the inner edge of the disk is set 
by dust sublimation and not the dynamical carving of 
a companion. The radius at which this occurs depends 
on the grain properties and the flux striking the surface, 
which in the case of LRLL 21 is dominated by the stellar 
flux (see discussion below) . Since the stellar flux is rela- 
tively constant on the timescales over which the infrared 
flux varies, the radius of the inner disk should be con- 
stant. The changing flux implies that the emitting area 
of the disk is varying, and if the radius of the inner disk 
is constant then the scale height must be varying. 

Since the 2-5/im excess is reprocessed stellar flux, the 
ratio of the excess to the stellar flux represents the cover- 
ing fraction of the inner disk, which we can use to quan- 
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Fig. 10. — Dereddened spectra of LRLL 21, along with the veiUng 
measurements (orange squares), 3.6 and 4.5fim photometry where 
available (blue triangles) and a stellar photosphere for comparison 
(blue line). Red dashed line is the best fit photosphere+blackbody 
emission for those days where we have spectra beyond 2.5^m. 
LRLL 21 shows a large range in the strength of the excess through- 
out our observations. 

tify the strength of the excess on different epochs. We 
calculate the total dust emission by integrating a black- 
body at the measured strength and temperature. For 
those days in 2009 when we do not have 2-5/im spectra 
we instead use the 3.6 and 4.5/im photometry to estimate 
the strength of the excess and assume a temperature of 
TdMst=1900K. The results are reported in Table [Hj for 
LRLL 21 the covering fraction varies from ~ 4% up to 
13%, while typical T Taur i stars have inner disks with 
covermg fractions of 12% (|Dullemond et al.ll2001[ ). On 
average the covering fraction of the inner rim of LRLL 21 
is lower than a typical T Tauri star, which may be due to 
advanced grain growth and settling at the very inner edge 
of the disk, although the strong silicate emission feature 
is inconsistent with significant grain growth where this 
flux arises. Another possibility is that the inner rim does 
not azimuthally extend all the way around the star, as 
was suggested earlier to explain the high measured dust 
temperature. If instead the inner disk has patches of op- 
tically thin material interspersed between optically thick 
dust, then the total covering fraction could be reduced. 
This may explain the periods when the emission is espe- 
cially low, which would otherwise require an extremely 
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Fig. 11. — Line measurements from 2008 (top two panels) and 
2009 (bottom two panels). Both the Pa/3 Une (A = 1.282^m) and 
the Br7 line (A = 2.166/im) are shown for each day. The solid 
line is the observed spectrum while the dashed line has had the 
photospheric absorption subtracted from the line. At the top of 
each panel is the KO WTTS standard. There are large variations 
in the emission lines, mainly on timescales of weeks. 

flat disk seen close to edge on. The lowest measured cov- 
ering fraction would correspond to only one-quarter of 
the inner disk being occupied by optically thick mate- 
rial, while the rest was either optically thin or devoid of 
dust. 

2.4.3. Gas Properties 

Based on eight total epochs of near infrared spectra, 
including four in 2008 over the course of 8 days and 
four in 2009 over the course of one month, we find large 
fluctuations in the Pa/3 and Br7 lines, often changing 
from absorption to emission fFig fTTj) . The emission line 
strengths of Pa/3 and Br7 reach a peak of -1.25 and - 
1.11 A respectively, with accretion rates varying from un- 
detectable to 10"^Moyr-i (Table E]). We place upper 
limits of a few times lO~^M0yr~-'^ on the accretion rate 
based on the uncertainty in the spectra on those days 
when no emission was detected. There are no day-to- 
day changes seen in the 2008 spectra but there are large 
changes fr om one week t o th e next, as is also seen in the 
2009 data. iDahml (|2008[ ) and lLuhman et al.l (|2003D mea- 
sure the EW of the Ha line to be -4.48 and -4.7A respec- 
tively, while our high-resolution MMT spectrum (Fig [T2]) 
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Fig. 12. — Optical spectra of the Ha line for LRLL 21. Dashed 
line is the median velocity of the cluster, while the two tick marks 
are placed at the velocities of the two binary components. 

shows an EW of -7.1 A, all above the sp ectral type de- 
pendent Ha EW threshold for accretors ([White fc Basril 
[200l . and our low-resolution spectrum measures -1.59A, 
consistent with large fluctuations in the accretion rate. 
In our high-resolution Ha spectrum the full velocity 
width at 10% maximum height is 422 km /sec, above the 
270 k m/sec threshold for active accretors (jWhite &: Basril 
12003( 1 ■ ahhough this may be influenced by the possible 
presence of a com panion. The high-resolution profiles of 
the Ha line from iDahml (j2008[ ) and our MMT spectra 
show strong central absorption, possibly due to two sep- 
arate Ha emission lines from the two components of the 
binary system as well as possible wind absorption. The 
possible binary nature of this system prevents us from us- 
ing the resolved line width to der ive the ac cretion rate, as 
has been done in other studies (jNatta et al. 2004). Our 
low-resolution optical spectrum, taken in 2009 when the 
Pa/3 line was in absorption, shows especially weak Ha 
emission compared to what is typically seen around ac- 
tively accreting stars and is like ly entirely chromosph eric 
rather than accretion emission (I White fc Basril[200l . al- 
though it is possible that unresolved absorption leads to 
a very small EW even with ongoing accretion. This very 
weak accretion measure is similar to the behavior of the 
infrared lines, which show no emission in some epochs, 
possibly also due to unresolved absorption. Even during 
the most active accretion, the ratio of accretion to stellar 
luminosity (Lacc/L*=0.06) is small. 

2.5. LRLL 58 
2.5.1. Stellar Properties 

LRLL 58 is a Ml. 25 s tar with a lumi nosity of O.72L0 
and a radius of 2.IR0 (jLuhman et al.l[2003i) . Its mass 
is ~ 0.7Mo based on its position on the HR dia- 
gram relativ e to th e iSiess et all (120001) 3 Myr iso chrone. 
iCohen et all (|200l and iCieza fc Balibeii (pOOel ) detect 
periodic variability in the optical with a period of 7.4 
days, consistent with the rotation of star spots across 
the stellar surface. The I band flux is roughly constan t 
overyears (1=14 10.14.32 14.18 fromlCohen et al.l ()2004[ ) , 
ICieza fc Balibol (I2OO6I) . ILittlefair et al.l ([20051) respec- 



tively) suggesting that there are no large fluctuations in 
stellar flux. From our hectochelle data we measure a ra- 



dial velocity of 14.2 km/sec and an upper limit on the 
rotational velocity of 15.0 km/sec (Table ITO)) . Based on 
two epochs of near-infrared spectra, its extinction is con- 
stant at kv=iA (Table El). Unlike LRLL 21, there is no 
substantial variation in the stellar flux of LRLL 58. 

2.5.2. Infrared Variability 

In the IRS spectra, LRLL 58 displays no variations 
(< 5% for A < 14^m and < 20% for A > 14^m), 
although the long-wavelength end has large systematic 
uncertainties due to uncertain background subtraction 
(Figure [T5)) . These systematic uncertainties are the likely 
source of the mismatch between the SL and LL spectra 
at 14/im, which we correct for by decreasing the LL spec- 
tra by a constant multiplicative factor when comparing 
different epochs. The 24/zm photometry shows signifi- 
cant fluctuations (14%), which fall below our detection 
limit in the IRS spectra, on predominately weekly and 
longer timescales (Fig lTl)) . The five consecutive days of 3- 
8/im photometry show a steady increase in flux, while the 
3.6 and 4.5/im monitoring confirms that large (0.3 mag), 
daily fluctuations occur, with no evidence for periodic- 
ity. When comparing the IRS spectra with the 2-8/im 
photometry there appear to be very large fluctuations 
that have not been captured during any single monitor- 
ing campaign. The 3-8/im photometry shows a distinct 
excess in all but one epoch, while the IRS spectra is con- 
sistent with almost no excess at these wavelengths. Only 
at 8/im are the photometry and the spectra consistent 
with the same flux level. In the 3-8/Ltm monitoring there 
is a distinct wavelength dependence: as the wavelength 
increases the size of the fluctuation decreases. This is 
consistent with the SED variations between the IRS spec- 
tra and the photometry suggesting that the SED pivots 
at A = 8/im, similar to LRLL 21 and 31. Significant veil- 
ing is present in the 0.8-2.5/im spectra (r/^ =0.46,0.48, 
Table [TTlFig I15ll6p . As with LRLL 2 we cannot con- 
strain the exact source of this excess, but it is consistent 
with hot dust emission (> 1200K) from close to the star. 
As with LRLL 21 we can estimate the covering fraction 
of the hot dust using our near-simultaneous near-infrared 
spectra and Spitzer photometry on 9 Oct 2009. We find 
a covering fraction of ~13% (Tabic [TT|) . which is compa- 
rable to what is seen around typical T Tauri stars. This 
represents the inner disk at its maximum flux, while there 
are some epochs when the 2-8/im flux drops almost to the 
photosphere, suggesting a very small covering fraction. 

2.5.3. Gas Properties 

During our two near-infrared observations in 2009 sep- 
arated by roughly one month the Pa/3 and Br7 emission 
appears to change. FigllTlshows that in the second epoch 
the lines are slightly stronger, although the noisy contin- 
uum makes it difficult to calculate accurate EWs. An ad- 
ditional complication is that our Ml WTTS standard JH 
108, which is used to estimate the photospheric absorp- 
tion, appears to show emission near the Br7 line (Fig|17p. 
making it difficult to estimate the emission for these lines. 
The equivalent width of the Ha line varies from -9A 
(|Herbiglll998h to -20.8A from our low-resolution spectra 
and is -14.9A (full-width at 10% of 240km/sec) in our 
high-resolution spectrum (Fig ITSl) , indicating that there 
have been large fluctuations in the accretion flux. The 
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Fig. 13.— Same as Figurc[2]but for LRLL 58. 
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Fig. 14.— Same as Figure |3] but for LRLL 58. 

small Ha line strength measured bv lHerbi^ (|1998[ ) is close 
to the boundary between a ccreting: and non-accreting 
sources () White fc Basri|[2003l ). suggesting that the accre- 
tion rate can drop to very low levels. Our high-resolution 
spectrum displays a broad line width indicative of ongo- 
ing accretion. Based on our measured accretion rate from 
the infrared lines (M 5 x 10~^MQyr~^), the accretion 
luminosity is a small proportion (Lacc/L*< 0.1) of the 
total luminosity that is responsible for heating the disk. 
The accretion rate derived from the width of the Ha line 
(|Natta et al.ll2004i ) is M - 3 x 10-"MQyr-\ which con- 
firms that the accretion luminosity is very small. Given 
the variable accretion flow and the systematic uncertain- 
ties in these accretion measures, it is difHcult to evaluate 
the consistency in the accretion rate measured by the 
infrared lines and the resolved width of the Ha line. 

2.6. LRLL 67 
2.6.1. Stellar Properties 

LRLL 67 is a MO . 75 star with a luminosity of O.48L0, 
a radius of l.GR© (jLuhman et al.. .20033 . and a mass of 



~ 0.6Mq (jSiess et al.ll2000[ ). iDahml (|2008[ ) measured 
the radi a,l velocity to be 1 5. 09 km /sec, with vsini=13.59 
km/sec, iNordhagen et al.l ()2006D find Vr=16.9 km/sec 
with vsini<11.0 km/sec and our hectochelle data mea- 
sure Vr=15.0km/sec and vsini<15.0 km/sec f Table ITO)) . 
The lack of significant variation in the radial velocity, 
sugge sts that there is n o massive companion in this sys- 
tem. iCieza fc Balibei] ([2006) did not find evidence for 
periodic optical fluctuations. We do not have any near- 
infrared spectra of this star with which t o estimate the 
extinction and we adopt the value from iLuhman et al.l 
(|20M) (Ay=2.4). As with LRLL 2 and 58, there is no 
evidence for large fluctuations in the stellar flux. 

2.6.2. Infrared Variability 

In the IRS spectra, LRLL 67 displays no variations 
(< 3% at A < 14^m and < 10% at A > 14^m) (Fig [HI). 
The 24/im photometry display marginal fluctuations at 
the 2(7 level, which are consistent with a constant flux, 
although there is a deviation between the 24^m photom- 
etry points and the long- wavelength IRS spectra (Fig [20)) . 
The flve consecutive days of 3-8/zm photometry show no 
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Fig. 15. — Spcx spectra (black line) and photometry (blue points) 
for LRLL 58. Spectra have been scaled to match 2MASS photom- 
etry. A stellar photosphere (JH 108 + Kurucz model) shown in red 
has been reddened by the Ay measured on each night to LRLL 58 
and scaled to the short wavelength flux of the spectra. The mis- 
match between the spectra and the 2MASS photometry is likely a 
sign of variability between these epochs. 

fluctuations (< 5%), but the 3.6 and 4.5/im monitoring 
does show large (0.4 mag) variations with no periodicity. 
Most of the change during the monitoring occurs during 
a brightening event with a timescale of about a week. 
Prior to this the infrared flux varied by less than 0.2 
mag, which is still significant given the small uncertain- 
ties in these data. The measurements are consistent with 
larger fluctuations at short wavelengths, similar to LRLL 
21, possibly involving periods of quiescence. During the 
3.6,4. 5/xm monitoring the [3.6]- [4.5] color gets redder as 
the source gets brighter (Fig [T5| , in a way that is con- 
sistent with a changing flux level from warm dust close 
to the star, similar to LRLL 21. lEspaiUat et al. (2012) 
find that a model with no optically thick dust within 10 
AU is able to fit the infrared SED. The implications for 
the variability of the possible lack of optically thick dust 
close to the star on the variability are discussed in more 
detail in section 13.1.21 

2.6.3. Gas Properties 

We h ave n o near-infrared spectra of LRLL 67, but 
IDahml (pool measure an Ha equivalent width of - 
28.32 A with an accretion rate of 2 x lO~^M0yr~^ 
while iLuhman et al.l ()2003() measure an Ha line strength 
of -35A and we find EW=-31.4A (full-width at 10% 
of 310km/sec) in our high-resolution optical spectrum 
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Fig. 16. — Dereddened spectra of LRLL 58, along with the veiling 
measurements (orange squares), 3.6 and 4.5/xm photometry where 
available (orange triangles) and a stellar photosphere for compar- 
ison (blue line). LRLL 58 shows a strong excess throughout our 
spectra. 

(Fig [21]) . LRLL 67 appears to be actively accreting, with 
a variable accretion rate, whose accretion luminosity is 
small (Lacc/L*=0.06) compared to the stellar luminos- 
ity. The accretion rate derived from the width of the Ha 
line UNatta et al.ll2004[) is M ~ 1 x IQ-'^^M^yr-^ , which 
confirms that the accretion luminosity is very small. 

2.7. LRLL 1679 
2.7.1. Stellar Properties 

LRLL 1679 is the faintest star in our sample and was 
cla ssified as a M3 . 5 star with L=O.21L0 and R*=1.3R0 
bv iMuench et all (I2007D . W e estimate its mass as ^ 
O.3M0 using the iSiess et all ()2000) isochrones. Based 
on five epochs of near-infrared spectra (three in 2008, 
two in 2009) we find that its extinction is constant at 
Ay =5. 8 (Table [5]). There is no mention of periodic op- 
tical fluctuations in the literature, although this source 
may have been too faint for these surveys or outside the 
field of view (at Dec=31:58:255 it is south of the cluster 
center). 

2.7.2. Infrared Variability 

In the IRS spectra, LRLL 1679 displays large fluctu- 
ations (^ 50%) on timescales of months, but no weekly 
changes (< 5% for A < 14^m and < 20% for A > 14/.tm) 
(Fig [22]). The SED appears to pivot, similar to LRLL 2, 
21, 31 and 58, at A = 13/im. The 24/xm photometry dis- 
plays fluctuations (50%) that are consistent with those 
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Fig. 18. — Optical spectra of the Ha line for LRLL 58. Dashed 
line is the mean velocity of the cluster. 



Fig. 17. — Same as Figure [S] but for LRLL 58. There is some evidence for a change in the line strengths between the two epochs, and 
ongoing accretion in both epochs, although the uncertainties in the derived accretion rates are large. When subtracting by a photospheric 
standard we consider both a WTTS and a giant from the IRTF spectral library, since the WTTS seems to exhibit Br7 emission. 

much smaller than expected for a typical T Tauri star. 
The small emitting area of the inner disk may be due 
to either substantial grain growth and settling reducing 
the scale height of the inner disk, or non-axisymmetric 
distribution of optically thick material at the inner edge. 

2.7.3. Gas Properties 

During our five epochs of near- infrared observations, 
three in 2008 over the course of a week and two in 2009 
separated by a month, we find no evidence for Pa/3 or Br7 
emission (Fig [^5)1 . We can place upper limits on the ac- 
cretion rate of 10~^ MQyr~^ (Table[7|). Our low resolu- 
tion optical spectrum taken in 2009 during the 3.6,4.5/im 
monitoring measures an EW of the Ha line of -3.3A, 
which is below the accretion boundary for a star of this 
spectral type ( White & Basri 2003). The high- resolution 
Ha spectra has an EW of only -1.7A (full-width at 10% of 
120km/sec) and shows very weak, narrow emission con- 
sistent with completely photospheric emission (Fig [26l) . 
In no epoch do we detect evidence for ongoing accretion 
in this star. 

2.8. Summary 

Our sample represents a wide range of stellar masses, 
from O.3-3.5M0. Except for LRLL 21, none of the sample 
show evidence for large changes in the stellar flux, either 
due to occultation events or hot spots rotating across 
the surface of the star. In previous optical monitoring 
surveys of IC 348 some do show small fluctuations due 
to cold spots rotating across the surface of the star. In 
the case of LRLL 21, changes occur from one year to the 
next, but not on daily or weekly timescales. 

The exact size, wavelength dependence and timescale 
of the infrared fluctuations varies from source to source, 
but overall all six sources show variability at all infrared 
wavelengths. Some show very large flux changes, up to 
50%, while others show much smaller variations. Most 
show fluctuations in which the short wavelength flux de- 
creases when the long wavelength flux increases. Some 
sources (LRLL 2,31) show a correlation between the sil- 
icate flux and the long-wavelength flux, while another 
source (LRLL 1679) does not, which may reflect differ- 
ences in the geometry of the dust and the location of 



seen in the IRS spectra (Fig [23]). Despite the large fluc- 
tuations in the continuum, there is no change in the flux 
of the silicate feature (although it is the weakest in our 
sample). The five consecutive days of 3-8^m photome- 
try show no variation (< 5%), consistent with the longer 
variability timescale. This is also confirmed in the 3.6 
and 4.5/im monitoring, which shows a steady increase in 
the flux of 0.1 mag over the course of 40 days with no 
evidence for periodicity. The change in [3.6]-[4.5] color 
between the 3-8/im photometry and the 3.6,4.5/zm mon- 
itoring, where the source gets redder as it gets brighter, 
is consistent with optically thick dust in the inner disk 
getting brighter with time, although there is no evidence 
of veiling in the 0.8-2.5/zm spectra (Fig [24|l . The short- 
wavelength excess can be very small, with the measure- 
ments almost at the level of the photosphere in some 
epochs, and it is possible that when we obtained our near- 
infrared spectra the excess was below our detection limit. 
The non-zero [3.6]-[4.5] color suggests that there is some 
hot dust close to the star creating a weak excess. As with 
LRLL 21 and 58 we can estimate the covering fraction 
in 2009 based on our nearly simultaneous near-infrared 
spectra and Spitzer photometry. On Oct 9 and Nov 4 we 
find covering fractions of ^ 2 — 3% (Table fTTj) . which is 
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Fig. 20.— Same as Figuro|3]but for LRLL 67. 
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Fig. 21. — High-resolution optical spectrum of the Ha line for 
LRLL 67 normalized to the continuum. Dashed line is the mean 
velocity of the cluster. 



the silicate emission zone, in particular whether or not 
the silicate emission zone can be experience the same 
shadowing as the outer disk. The timescales range from 
weeks to months, with no evidence for long year-to-year 
changes, suggesting that we are not observing the rem- 
nants of a single burst or the slow removal/addition of 
dust to the disk. None of them show strong evidence for 
periodicity in the warm mission data, which is sensitive 
to P=2-25 days. There is some evidence for a depen- 
dence of the location of the pivot in the SED on mass 
with the wavelength of the pivot decreasing with increas- 
ing stellar mass, although it is more likely that this pivot 
wavelength depends more on the exact d ust geometry of 
the s ystem rather than the central mass (jEspaillat et al.l 

Eonl). 

In most of the sources there is evidence for variable 
accretion, which at some epochs dips below our detection 
limits. The only exception is LRLL 1679 in which there 
is no evidence of ongoing accretion above our detection 
limits at any epoch. Overall the accretion luminosity 
is small relative to the stellar luminosity even when the 
stars are actively accreting. 
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Fig. 22.— Same as Figurc[2]but for LRLL 1679. 
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Fig. 23.— Same as Figuro|3]but for LRLL 1679. 

3. VARIABILITY PHENOMENOLOGY 

Given all of this information, we seek to address a num- 
ber of questions: Is the variability in all of these systems 
caused by similar perturbations to the disk? Is the type 
of perturbation common among cluster members? 

3.1. Is the origin of infrared variability universal? 
3.1.1. Comparison with LRLL 31 

There is a bewildering array of infrared variability for 
our six sources. Some show large infrared fluctuations 
with a clear wavelength dependence in every epoch, some 
do not. Some (LRLL 58,67,1679) show almost no ex- 
cess emission for A < 8/im in some epochs, while others 
(LRLL 21) always have a significant excess. In Table [T2] 
we summarize the data presented in the previous section 
regarding the variability of the stellar flux, accretion flux 
and infrared flux, which we now use to determine if these 
sources can be included together as one coherent type of 
object. 

As a starting point we consider the mo del for the 
variability of LRLL 31 (|Flahertv fc MuzeroIIi 120101: 
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iFlahertv et al.l [20111 1 whi ch relies on a vertica l pertur- 
bation to the inner disk. lEspaillat et al.l (|2011l ) employ 
a similar geometry, with a much more detailed model 
of the dust structure and the radiative transfer, to suc- 
cessfully model the infrared fluctuations of a sample of 
nearby transition disks, suggesting that this model may 
be applicable to many sources. The key component is 
optically thick material close to the dust sublimation ra- 
dius that experiences some sort of vertical perturbation 
causing it to rise out of the midplane. The consequences 
of this model are: 

• An SED that varies at all infrared wavelengths and 
appears to pivot around a certain wavelength be- 
tween different epochs. As a vertical perturbation 
rises out of the inner disk it produces more short 
wavelength flux while also shadowing the outer 
disk, causing the long-wavelength flux to decrease. 
In the case of the LRLL 31, the flux shortward of 
8/^tm would decrease when the long-wavelength flux 
increased. The pivot point may change depending 
on the geometry of the disk (Espaillat et al...2011.) , 
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Fig. 24. — Spex spectra (black line) and photometry (blue points) 
for LRLL 1679. Spectra have been scaled to match 2MASS pho- 
tometry. A stellar photosphere (LkCa 21 + Kurucz model) shown 
in red has been reddened by the Ay measured on each night to 
LRLL 1679 and scaled to the short wavelength flux of the spectra. 
Very little excess is seen throughout the observations. Spectra have 
been scaled by a constant factor for clarity. 
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Fig. 25. — Same as FigurelTTI but for LRLL 1679. We do not sub- 
tract the photospheric absorption since it is very weak and would 
only add to the uncertainty in the spectra. There is no evidence 
for either Pa/3 or Br7 emission in any epoch. 

but it is difficult to produce this wavelength de- 
pendence without a variable inner disk scale height 
(jFlahertv fc MuzeroUd |2010[ ). This perturbation 
may or may not be axisymmetric about the star. 

• The infrared flux changes on the dynamical 
timescale of the inner disk (typically days to 
weeks). These rapid variations can even be seen 
at the longest wavelengths, which is dominated by 
emission from material at a few AU from the star 
where the dynamical timescales are closer to years 
and the disk structure cannot substantially change 
in a matter of days. The long wavelength fluc- 
tuations are due to variable shadowing instead of 
in-situ structural changes. 

• The [3.6]- [4.5] color gets redder as it gets brighter. 
A change in the disk emission, rather than simply 
the stellar flux, can produce this effect since the 
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Fig. 26. — Optical spectra of the Ha hne for LRLL 1679. Dashed 
line is the mean velocity of the cluster. 

relative contribution of the disk emission, specifi- 
cally the inner disk which dominates at these wave- 
lengths, is larger at 4.5/^m than at 3.6/im. 

Overall each source displays some of the signs of a vari- 
able inner disk suggesting that all of these sources are 
experiencing the same type of variability. LRLL 2 dis- 
plays a SED that pivots around A = 6/im. LRLL 21 
has a light curve that is almost identical to LRLL 31; its 
[3.6]- [4.5] color gets redder as it gets brighter and direct 
observations of the emission from the inner disk show 
that it is optically thick and the size of the emitting re- 
gion gets larger while its temperature stays the same. 
The 24/xm photometry shows that LRLL 21 is variable 
at this wavelength. LRLL 58 shows a pivot in the SED 
at A = 8/xm when comparing the photometry and the 
IRS spectra, and the 24/im photometry shows significant 
fluctuations. It also shows variations on short timescales, 
consistent with the dynamical timescale of the inner disk 
(5 days). LRLL 67 shows rapid fluctuations consistent 
with the dynamical timescale of the inner disk (4 days) 
and the [3.6]- [4.5] color gets redder as it gets brighter. 
LRLL 1679 clearly displays a pivot in its IRS spectra at 
A = 13/im and there is some evidence that the [3.6]- [4.5] 
color gets redder as it gets brighter. The small change 
in optical flux on short timescales in all of these sources 
suggests that the irradiation of the disk does not change 
appreciably. If the inner disk radius is set by dust sub- 
limation and thus proportional to the stellar and accre- 
tion luminosity, it is not changing in these objects. The 
rapid infrared fluctuations imply changes in the emitting 
area of the disk, but since the radius of this dust does 
not change then the scale height must be variable. That 
none of them show all of the hallmarks of a variable inner 
disk may be a sign of differences in the exact geometry of 
the disk (amount of optically thin dust in the inner disk, 
location of the outer disk, grain size/composition) rather 
than differences in the perturbation. In general, the sim- 
ilarities in behavior suggest roughly similar underlying 
causes for the variability with the possible exception of 
LRLL 67. 

3.1.2. The case of LRLL 67 

One key tenet of this model for the infrared variability 
is an optically thick inner disk and in most of the tran- 
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sition disks we have strong evidence for such material. 
LRLL 21 and 31 show direct evidence for optically thick 
inner disks based on the shape of the 2-5/im excess. The 
strength of the veiling observed at K band in LRLL 58 
is also consistent with strong emission for an optically 
thick inner disk, rather than the weak emission from an 
optically thin disk. LRLL 2 and 1679 show a pivot when 
comparing IRS spectra, while LRLL 58 shows a pivot be- 
tween the spectra and photometry, which is difhcult to 
explain without shadowing from an optically thick inner 
disk. For LRLL 67 the evidence for an optically thick 
inner disk is weaker than in any other star. Its SED is 
reminiscent of sources such as GM Aur, DM Tau, CS 
Cha and Coku Tau/4, which have been modeled using 
only optically thin dust in t he inner d isk fEsp aillat et al.l 
[^07; D'Alessio et al. 2005; Calvet et al. 2005). We have 
no near-infrared spectra to look for a sign of hot dust, 
but the weak 3-8/im excess at all epoch s is difficult to 
explai n with optically thick dust. In fact. lEspaillat et"al] 
are able to fit the infrared SED without the use of 
any optically thick dust within 10 AU. If the inner disk 
is optically thin then changing the scale height will have 
no effect on the short wavelength flux, since the emission 
is proportional to the total column density of material 
and not its exact geometry. The inner disk will also not 
be able to shadow the outer disk, leading to no variation 
in the long-wavelength flux. 

LRLL 67 has a substantial accretion rate, indicating 
that some material must be traveling from the massive 
outer disk and onto the star implying that the inner disk 
is not completely devoid of material. The measured ac- 
cretion rate of 2 x 1O~^M0 yr~^ corresponds to a gas 
surface density of ~ 10 g cm~^ at 1 AU as suming a stan- 
dard a disk model (jD'Alessio et al.l[l999l ). If the gas to 
dust mass ratio is 100, then the d ust emission should 
be optically thick. iRice et all (|2006( ) find that the outer 
edge of the dust depleted gap can block some of the in- 
ward flowing dust resulting in a very high gas to dust 
ratio and a much lower dust mass in the inner disk, pos- 
sibly accounting for the weak dust emission despite the 
high gas density. Another explanation for the weak ex- 
cess is if the inner disk were radially optically thick while 
vertically optically thin. In this case it would be able to 
shadow the outer disk, while still producing very weak 
emission. This model has difficulty explaining the vari- 
ability at 3.6 and 4.5/^m with a changing geometry, al- 
though if the inner disk were optically thin then a change 
in its temperature could lead to substantial changes in 
the infrared flux. This variable heating coul d be due to 
an un seen companion on an eccentric orbit (|Nagel et all 
I2010D or chan ges in the accretion rate, although for LRLL 
67 the heating of the disk is dominated by the stellar flux. 

A completely optically thin inner disk would produce 
no shadowing o f the outer disk , and no 'seesaw' behav- 
ior in the SED. lEspaillat et"all ()20lH) find that the pre- 
transition disks (i.e. those that have optically thick in- 
ner disks) in their sample show the 'seesaw' behavior 
in the mid-infrared spectra, while the sources without 
optically thick dust in the inner disk do not. Of our 
evolved disks, LRLL 67 shows the least amount of evi- 
dence for a pivot in the SED or any 'seesaw' behavior 
in the IRS spectra, leaving open the possibility that the 
outer disk is not shadowed, although the difference be- 
tween the 24/Ltm emission and the IRS spectra is difficult 



to explain without shadowing of the outer disk. This long 
wavelength variability suggests that some fraction of the 
inner disk should be optically thick. Observations of the 
near-infrared excess to determine if it is optically thick, 
as has been done for LRLL 21 and 31, as well as detailed 
modeling of the SED are needed to better understand 
this source. 

While the other evolved disks all have strong evidence 
for optically thick material in the inner disk, they also 
all have epochs when the short- wavelength excess is very 
weak, which is difficult to explain with an optically thick 
disk. For LRLL 21 and 1679 we estimate covering frac- 
tions of 0.02-0.04 when the excess is at its weakest (Ta- 
ble[TT|). A perfectly flat disk, which is needed to minimize 
the flux, would need to be very close to edge on to pro- 
duce an excess this small. It is also possible that the 
inner disk is non-axisymmetric, with the majority of it 
optically thin with only thin trails of optically thick dust. 
The small covering fr action of these trails would cre - 
ate a very weak excess UDodson-Robinson fc Saiv3[20Tll) . 
LRLL 67 could represent an extreme case with only very 
thin trails of optically thick material, while the other 
stars have larger patches of optically thick material. In 
addition to vertical perturbations of the inner disk, there 
may also be variations in the azimuthal coverage of the 
optically thick material. In some epochs the azimuthal 
coverage could be very small, producing a minuscule ex- 
cess at short wavelengths. Further observations and SED 
modeling of these systems are needed to fully constrain 
the possibly non-axisymmetric distribution of material 
in the inner disk. 

3.2. Comparing Transition Disk Variables to Other 
Variables 

Do these stars represent a unique class of objects, or is 
the physical mechanism responsible for their variability 
common to other young stellar objects? We can start to 
answer this question by comparing the dust grain prop- 
erties in our evolved disks, as traced by the 10/^m silicate 
emission, to those in a less evolved disk population. The 
emission can be characterized by the continuum normal- 
ized peak flux (S(peak)) and the ratio of the continuum 
norm alized flux on two si d es of the feature (S'li.s/S'g.s) 
as m iKessler-Silacci et al.l (|2006f ). The wavelengths 11.3 
and 9.8/im correspond to the peak of the emission proflles 
from crystalline and amorphous silicates, allowing us to 
distinguish processing of the dust grains. In Figure [27l we 
compare the strength and shape of the silicate feature for 
each observation of the evolved disks with the t rend pre- 
viously observed bv IKessler-Silacci et al.l (|2006D and flnd 
that the dust properties of the evolved disks studied here 
are similar to those of other star s. A similar conclusi on is 
found when using t he method of Ma noi et al.l (j2011| ) and 
iFurlan et al.l ([20091 ) , which apply slightly different indices 
to trace the strength and shape of the feature, and com- 
paring to the Taurus, Chameleon and Ophiuchus cluster 
members studied in those papers. Although our sample 
consists entirely of evolved disks, and by dcflnition the 
structure of the dust within the disk is different from a 
typical T Tauri star, there is no clear distinction in their 
dust properties. The previously mentioned studies have 
also found that dust properties do not vary strongly be- 
tween normal and evolved disks. 

The infrared fluctuations of these evolved disks origi- 
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Fig. 27. — Strength versus shape of the 10/jm silicate feature for 
the evolved disks. Definitions of the axes are from Kessler-Silacci 
et al. 2006. The dashed Une shows a fit to stars from c2d and 
demonstrates the correlation between strength and shape found 
for typical stars. The evolved disks tend to lie slightly above this 
relation, but are not substantially different from other T Tauri 
stars. 

nally drew us to this sample, and they can be used as a 
point of comparison with other young stellar objects. We 
have used our IRAC and MIPS photometry to compare 
the size and timescale of the infrared variability between 
the transition disks and the rest of the young stars in IC 
348. One way to trace the typical size of the infrared 
fluctuations i s with the rms of the p hotometry. The rms 
is defined as (jCarpenter et al.ir2001| ): 



[rmsY 



(n 



(1) 



where n is the number of observations, rrii is the magni- 
tude on a given epoch, Wi = l/crf and Oi is the photomet- 
ric noise. In Figure US] we show the rms of the fluctua- 
tions at every wavelength for which we have photometry, 
with the size and color of each point scaled to the slope, 
a, of the IRAC SED, defined as \F\ oc A", since the fre- 
quency of infrared variability a ppears to depend on the 
strength of the infrared excess (|Morales-Calder6n et al.l 
120111 ). For the [3.6] and [4.5] photometry we only include 
the warm mission photometry since any systematic un- 
certainties in the photometry may differ between the cold 
mission and warm mission. For the [5.8] and [8.0] pho- 
tometry we exclude the c2d photometry because the dif- 
ferent observing strategy prevented them from measuring 
accurate photometry of many of the brightest members. 
We also mark the transition disks in IC 348 identified by 
[MuzeroUc ct al. (2010) for a comparison between evolved 
disks. The rms for the evolved disks studied here does 
not stand out from the rest of the cluster members at any 
wavelength. The number of epochs of data available at 
each wavelength varies (38 for [3.6] and [4.5], 6 for [5.8] 
and [8.0] and 10 for [24]) but the conclusion is the same. 
The fluctuations of our sample are larger than average, 
but this is mainly because these sources were chosen as 
the evolved disks that showed the largest fluctuations. 

The similarly sized fluctuations suggests that struc- 
tural changes in disks, similar to those seen in the 
evolved disks studied here, are common among young 
steUar objects. Roughly 60% of the IC 348 cluster mem- 



bers with an infrared excess show variability at 3.6 and 
4.5jHni larger than 0.04 mag (Flaherty et al in prep.). 
iMorales-Calderon et al] ()201lD find that 70% of the stars 
in Orion with an infrared excess are variable at 3.6 and 
4.5/im with amplitudes similar to those studied here. 
The amplitude of the fiux variations is very similar be- 
tween the evolved disks studied here and the other vari- 
able disks, which make up a majority of the members 
with disks. Despite the similarities in rms amplitude it is 
still possible that the variability mechanisms differ. Fur- 
ther observational evidence is needed to probe whether 
the normal disks are experiencing vertical perturbations 
of the inner disk, rather than other structural changes. 

There is evidence for varying, non-axisymmetric struc- 
ture in inner disk s around other stars in the literature. 
iHutchinson et all (|l994) find a wavelength dependence 
in the variability of UX Ori and AK Sco that can be ex- 
plained by changing the sca le height of the inner disk 
^Flahertv fc Muzerolld l20Tol ). lJuhasz et all (|2007l) see 
infrared variability in SV Cep that they are able to 
model successfully with a disk whose inner scale height is 
rapidly varying. None of these objects are characterized 
by transition disk SEDs. 

Other systems also show evidence for non- 
axisymmetric structure in the inner disk, based on 
the occasional obscurat i on of the star by this material. 
IMorales-Calderon et al.l (|2011[ ) find a significant fraction 
of the variables in Orion display rapid dips in their 
light curves due to brief stellar occultations by the 
disk. The COROT satellite monitored the young 
cluster NGC 2264 for 23 straight days and found that 
28% of the young s tellar objects exhibit this behavior 
(|Alencar 'et^mm . AA Tau is a well-studied example 
of a star with periodic occultations due a warp in the 
inner disk (Bouvicr ct al. 2007a). The structure that 
causes these occultations may be long-lived, or it may 
rapidly appear and disappear on an orbital timescale. In 
fact UX Ori, which shows strong evidence for an inner 
disk with a variable scale height based on its infrared 
fiuctuations, is the canonical member of the UX-Ori 
class of variables t hat are characterized by brief drops in 
their optical flux (jHerbst fc: Shevchenkol il999). similar 
to those seen in Orion, NGC 2264 and AA Tau. For 
this class of objects the disturbances in the inner disk 
are most likely non-axisymmetric in order to produce 
periodic occultations. The frequency of observations 
that are consistent with variable non-axisymmetric inner 
disk structure suggest that this is a fundamental feature 
of protoplanetary disks. 

4. WHAT COULD BE CAUSING THESE VARIATIONS? 

Our multi-wavelength data tracing more than just the 
dust emission allows us to put some constraints on the 
physical cause of the infrared variability. We use these 
data to eliminate accretion and winds, but companions 
or magnetic flelds are still plausible. 

4.1. Variable Irradiation of the Disk 

One possible explanation for the observed infrared 
variability is that the illumination of the disk changes 
rapidly. As the flux striking the disk increases the tem- 
perature of the disk rises. If the inner edge is set by dust 
sublimation then the location of the inner edge will move 
outward, increasing its emitting area and hence flux. The 
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increased heating may also raise the temperature of the 
midplane resulting in an increase in the disk scale height, 
shadowing the outer disk. The illumination of the disk is 
dominated by the stellar and accretion luminosity, which 
we can trace using our near-infrared photometry and our 
measurements of the Pa/3 and Br7 lines as well as with 
measurements from the literature. The observations of 
the stellar and accretion flux for these stars are not as 
dense as the infrared coverage, but we can still look 
for variability in these properties. Our contemporane- 
ous near-infrared spectra, along with near-infrared and 
mid-infrared photometry will allow us to look for corre- 
lations among the stellar flux, accretion flux and infrared 
flux although it may be challenging to establish a causal 
relationship. 

In terms of stellar flux, none of the six evolved disks 
show evidence for strong fluctuations on timescales simi- 
lar to that of the infrared variability. LRLL 2, 58, 67 and 
1679 display, at the most, small changes due to cold spots 
rotating across the face of the star. These fluctuations 
are not large enough to substantially change the temper- 
ature of the disk. LRLL 21 and 31 show evidence for 
year-to-year changes in stellar flux, which is much longer 
than the timescale for the infrared variability. For both 
stars we have contemporaneous J band and 3.6,4.5/im 
photometry in 2009 and we find that the J band flux 
stays relatively constant despite the large changes in the 
infrared flux. 

The accretion rate does show large variations in most 
of these stars, with the only exceptions being LRLL 2, 
which has large uncertainties on its accretion rate, and 
LRLL 1679, for which we do not detect any strong on- 
going accretion (< IG^^M© yr~^). Despite the accre- 
tion rate variability, the average accretion luminosity is 
small compared to the stellar flux and is not an im- 
portant source of heating for the disk surface. In none 
of the evolved disks is the accretion luminosity larger 
than 10% of the total flux striking the disk. In order 
to change the scale height by a factor of flve, the irradi- 
ating luminosity would hav e to change by a factor of 9 
(jPuUemond fc MonieiH[2010l ). which is much larger than 
the observed variability. This indicates that even though 
the accretion luminosity changes significantly, the heat- 
ing of the disk, and hence its scale height, will not vary 
at an observable level. 

For LRLL 21 and 31 we have enough contemporaneous 
data in 2009 to look for any correlation between the ac- 
cretion rate and infrared flux. Both of these stars show 
similar infrared light curves, with a sharp drop during 
the first week of observations followed by a slow increase 
over the next four weeks. We were able to obtain Pa/3 
and Br7 measurements during the IR minimum as well 
as a few weeks later after the infrared fiux had increased 
for both stars. In LRLL 31, the Pa/3 line shows almost 
no emission during the infrared minimum while it has 
strong emission after the infrared flux has increased, sug- 
gesting a strong correlation between the gas and dust 
behavior. In LRLL 21, the Pa/3 line shows absorption 
during the infrared minimum (Oct 8) as well as 3 weeks 
later (Oct 31) after the infrared flux has increased by 0.7 
mag. One week later (Nov 4) the Pa/3 line is in emission 
even though the infrared flux has only increased by 0.1 
mag. It is possible that unresolved absorption on Oct 31 
makes the line appear very weak. The Hel line, which 



traces the inflow and outflow of gas near the star, shows 
very strong red-shifted absorption due to the accretion 
in this epoch, but not during any other epoch (Fig [29)) . 
There is also evidence for red-shifted absorption in the 
Pa7 line. This suggests that there may be additional 
absorption that suppresses the Pa/3 emission making it 
appear very weak. If this is the case then the Pa/3 flux 
may be correlated with the infrared flux as in LRLL 31. 
Even if the accretion rate and infrared flux are correlated 
it is unlikely that the accretion variability is the source 
of the infrared fluctuations. As discussed earlier, the ac- 
cretion luminosity is small compared to the stellar flux 
heating the disk. It is also unlikely that a variable ac- 
cretion flow through the disk leads to the variable disk 
emission because the timescale for an over-dense clump 
of dust and gas to move through the disk is the viscous 
timescale, which is on the order of hundreds of years at 
the location of the inner disk edge for the objects. 

4.2. Winds 

A highly variable wind, carrying both gas and dust 
out o f the disk, could lead to variable infrared emis- 
sion ()Konigl fc Pudritd [2000t Mnkovfc & Jurkfc 2007). 
As material is lifted out of the disk its emitting area in- 
creases and it may shadow the outer disk, creating the 
observed infrared variability. It is also possible that the 
wind is asymmetric about the mi dplane, likely du e to a 
misaligned stellar magnetic field (|Lovelace et al.l 12010). 
If the mass flow in the wind is larger on one side of the 
disk than the other, then the asymmetric loss of angu- 
lar momentum will lead to warping of the inner disk. 
Both of these models require the wind to be launched 
outside of the dust su blimation radius, which is more 
typical o f a disk wind (|Konig 1 fc Pudritd I2OOOI ) than an 
X-wind (|Shu et al.l 119941 ). As discussed below, the co- 
rotation radius, which is approximately the radius at 
which an X-wind is launched, lies within the sublima- 
tion radius for three of our evolved disks (LRLL 2, 21 
and 31) preventing this type of wind from directly af- 
fecting the dust. We can trace the wind using the 
blue-shifted absorptio n / emission of the Hel (1.08/im) line 
(jEdwards et al.l 12006( 1 . When a clump o f material is 
lifted out of the disk midplane by a wind ([Suzuki et al.l 
|2010[) the clump needs to be optically thick to stellar 
photons to shadow the outer disk and produce the ob- 
served wavelength dependence of the variability. This 
sets a lower limit on the column density of material in the 
clump of Ay ^ 1, which corresponds to log(/V^f) ^ 21 
(|Drainel I2003D . A wind with this column density that 
arose from ~0.5AU and was able to remove material in 
a week would have an outflow rate of ~ 1O~^M0 
yr~^. This is comparable to the accretion rates detected 
in the evo lved disks, and should produce detectable Hel 
emission ()Kwan fc Fischer] 120 lit ) . In LRLL 2, 58 and 
1679 we see no evidence for absorption or emission in 
the Hel line (Fig^. Only LRLL 21 and 31 show blue- 
shifted emission in some epochs, but not nearly at the 
strength relative to the continuum expected for such a 
high column density of material. This suggests that any 
wind that is present in the system is not dense enough 
to shadow the outer disk. 

An asymmetric wind does not require the density in 
the wind to be as high, since it is the warped disk that 
is doing the shadowing of the outer disk, rather than 
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the wind itself. With this model we would expect the 
strongest blue-shifted absorption/emission, correspond- 
ing to the strongest wind, to occur during the epochs 
when the infrared flux, and hence the warping of the 
disk, is the largest. For LRLL 21 in 2009, when we have 
contemporaneous measurements of Hel and infrared flux, 
there is no sign of a correlation between Hel and the in- 
frared flux (Fig EH). In LRLL 31, in 2009, the strongest 
blue-shifted absorption occurs during the minimum of 
the infrared light c urve, which is the opp osite of the ex- 
pected correlation ()Flahertv et al.|[2011|) . This suggests 
that an asymmetric wind is not the source of the observed 
infrared variability. 

4.3. Companion 

A companion on an orbit that is misaligned with 
the disk will lead to a warp whose height varies as 
the companion drags m aterial out of the midplane 
(jFragner fc Nelsonl l2010t ) . If a companion were copla- 
nar with the di sk, then it could lead to variable heat- 
ing of the disk ()Nagel et alll2010[) or variable mass flow 
through the disk in the ga p surrounding the planet 
(jArtymowicz fc Lubowl[T996l ) . The fluctuating mass flow 
through a gap could lead to variable shock heating of 
the inner disk which would then change the scale height 
of the inner disk. The perturbation from such an ob- 
ject would occur on a timescale equal to the compan- 
ion's orbital period, which could be as little as a few 
days. There is evidence for radial velocity variations in 
the optical spectra of LRLL 21 and 31, which are con- 
sistent with the presence of a massive companion. We 
can put some constraints on the possible location of a 
companion based on two observations: (1) the existence 
of hot dust, responsible for the 2-5/im excess and (2) 
the lack of periodicity in the infrared for P=2-25 days. 
The flrst constraint prevents a companion from being lo- 
cated within 0.3a < r < 1.7a where a is the location of 
the inner dust, or else the in ner disk would be removed 
(jArtymowicz fc Lubowll 19961 ). It is possible the the inner 
disk has a finite radial extent, in which case the compan- 
ion is excluded from 0.3ao < r < 1.7ai where oq and ai 



are the inner and outer edges of the inner disk. We as- 
sume that the inner disk is infinitesimally thin in order to 
place the most conservative constraints on the location 
of a companion. We estimate the location of the dust 
destruc tion radii assuming Ti„,n,p.r= 1500 K and ISM like 
grains pullemond fc Monieiil [20101 ) . Here we have as- 
sumed that every system has optically thick dust at the 
dust sublimation radius. For LRLL 21 and 31 we have 
direct measurements of this excess, while in LRLL 2, 58 
and 1679 we have indirect evidence based on the strength 
of the K-band veiling and the 'seesaw' behavior in the 
IRS spectra. As discussed earlier, the evidence for opti- 
cally thick dust close to LRLL 67 is weak, although we 
include it here for completeness. We use ISM-like grains, 
but note that a change in the grain sizes and composi- 
tion can lead to substantially different estimates of the 
dust sublimation radius Espaillat ct al. (2012). The large 
spectral type range for our sample results in a large range 
of inner disk radii (Table [HI), from 0.76AU for LRLL 2 
down to 0.03 AU for LRLL 1679. The second constraint 
excludes companions from orbits with periods between 2 
and 25 days, assuming that a companion would perturb 
the material on every orbit. In Figure [30[ we indicate the 
radii for which a companion is excluded in these evolved 
disks. For LRLL 2 we can place the tightest constraints, 
and exclude a companion from 0.05-1.3AU. 

Two important, and related, caveats to the above anal- 
ysis: (1) dust can extend closer to the companion if it 
is in streamers rather than a uniform wall and (2) the 
effects of a companion on the disk may not be peri- 
odic. Simulations of the clearing of a disk by a com- 
panion find that small streamers of material can ex- 
tend from the undistur bed outer disk to the companion 
(jArtvmowicz fc Lubowl[T9 96!) . The limits used above for 
the location of the disk refer to the undisturbed outer 
disk and streamers can fill in some of the space that has 
been mostly cleared out by the gravitational influence 
of the companion. If these streamers are optically thick 
then they could create the small excess that is observed 
in some of our sources, due to the very small covering 
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fraction of this material. The niodels also show that 
the structure of these streamers may vary stochastically, 
with no strongly periodic component. This runs contrary 
to our assumption that a companion would perturb the 
disk on every orbit. The above analysis is still applicable 
to a companion on an inclined orbit, but breaks down 
when dealing with optically thick streamers of material 
that can arise from a companion being embedded within 
the disk. More direct tracers are needed to rule out or 
confirm the perturbation of the disk by a companion. 

4.4. Magnetic Fields 

Magnetic fields play an important ro le in channel- 
ing ma terial from the disk onto the star ([Bouvier et al.l 
l2007bl) and possibly for driving the turbulence within 
the disk that leads to the inward flow of material 
(jBalbus fc Hawle"vlll991| ). Many theoretical models pre- 
(lict that the structure of the magnetic field in these 
young systems should not be constant and these fluc- 
tuations may be the source of the vertical perturbations 
of the inner disk. Magnetic fields in the turbulent disk 
may occasionally become buoyant, lifti ng dust and gas 
off the disk on very sh ort timescales (Tu rner et al.l[2010l : 
IHirose fc Tilmeil [20TT1 ) . The stellar magnetic field, and 
where it intersects the disk, may also be important for 
variability. As material from the disk flows inward it 
eventually is loaded onto the stellar magnetic field lines, 
is lifted out of the midplane, and free-falls onto the star. 
A highly dynamic interface between the stellar magnetic 
field and the disk may lead to rapid changes in disk struc- 
ture. The stellar magnetic field may expand and contract 
on weekly timescales, lead ing to a warping of the disk 
as the field expands tGoo dson fc Wingled 119991 ) . If the 
stellar magnetic field is mis aligned w ith the disk, then 
this could also cause a warp ()Lai fc Zh ang 2008). Simu- 
lations also show that the infio w of material aloiig fiel d 
lines is unstable and very patchy ()Romanova et al.|[2009[ ). 
If dust is entrained in this patchy inflow then it could oc- 
casionally create a large short-wavelength excess as it is 
lifted out of the midplane as well as shadow the outer 



disk. 

These last two models (tilted magnetic field and patchy 
infiow) require dust to extend inwards toward the point 
at which material is loaded onto the magnetic field lines. 
This radius is typically taken as near the co-rotation ra- 
dius, where the periods of the stellar and disk rotation 
are equal, since outside of this radii the angular momen- 
tum of the gas and dust is large enough that it would 
be flung out of the system rather than accreted onto the 
star. The innermost radius of the dust is set by the point 
at which the dust is hot enough to sublimate, which de- 
pends on the lumi nosity of the cen tral star, as well as the 
dust composition (jPullemond fc Monier„201QV Table [HI 
lists estimates of the co-rotation radii, calculated using 
either measured rotation periods or the observed range 
of rotation periods for pre-main sequence stars in IC 348 
(jCieza fc Balibeiil2006D . as well as the dust sublimation 
radii, calculated earlier. For the low-luminosity sources 
(LRLL 58, 67, 1679) dust can extend close enough to 
the star to be loaded onto the stellar magnetic field lines 
and lifted out of the midplane. For the high-luminosity 
sources the dust disk is truncated at much larger radii, 
even if we assume that the disk is entirely composed of 
very large (> lOfim) grains, which can extend closer to 
the star than small (0.1/im) grains. We include estimates 
of dust sublimation radius using very large grains in order 
to put a lower limit on its location, despite the fact that 
the large silicate emission features in many of our IRS 
spectra are inconsistent with an inner disk dominated by 
large grains. 

Lifting dust out of the midplane would occasionally 
lead to higher extinction and decreased stellar flux co- 
incident with changes in the infrared flux as our line of 
sight passes through one of these dense clumps. This 
would be more prevalent the closer the dust approached 
the star. In none of our stars do we see evidence for an 
increase in extinction or a decrease in stellar flux associ- 
ated with the infrared flux variations. It is unlikely that 
for all of these six systems we have been lucky enough 
that our line of sight never intersects the dusty clump. 
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although the paucity of measurements of the extinction 
of stellar flux makes it difficult for us to say that the 
stars are never occulted. For LRLL 21 and 31, where we 
have the densest coverage of the extinction contempora- 
neous with the 2009 infrared monitoring, we do not see 
any correlation between these properties. This suggests 
that the coverage of the sky by the clump as seen from 
the star is small. Additional monitoring for extinction 
events associated with material being lifted out of the 
midplane would help to further constrain this theory. 

4.5. X-ray Flares 

Young stellar objects are known to exhibit stronger av- 
erage X-ray emission as well as stronger X-ray flare s than 
a ma in sequence star of similar mass (Fcigclso n et al.l 
120071 ) . X-ray flares last for up to 2 days and can substan- 
tially heat the disk (jWolk et al.l[2005l) . Recentlv iKe et al.l 
(|2012l) showed that flares similar in strength to the most 
powerful outbursts seen in Orion may even be respon- 
sible for the change in inner disk scale height. A large 
X-ray flare will dramatically increase the temperature 
and ionization of the dust and gas. The highly ionized 
dust will then be accelerated either inward or outward 
by the stellar magnetic field and this acceler ation will 
lead t o a change in the scale height of the disk. iKe et all 
(|2012D show that given a typical X-ray flare, amid as- 
sumptions about the disk structure, they can reproduce 
the 'seesaw' behavior on timescales of days to weeks that 
we observe here. Since the response of the disk to the X- 
ray flare is nearly instantaneous and lasts for days, con- 
temporaneous monitoring of the X-ray and infrared flux 
should be able to constrain this theory. And because the 
strength of the stellar magnetic field drops strongly with 
radius, those sources with dust close to the co-rotation 
radius should exhibit the largest fluctuations, assuming 
that different sources have similar X-ray flares. In our 
small sample we have objects that range from having 
dust very close to the co-rotation radius to objects with 
dust very far from the co-rotation radius and do not find 
any clear trend in the strength of the variability between 
these sources. 

5. CONCLUSION 

We have studied the mid-infrared variability of six 
evolved disk systems in the young stellar cluster IC 348. 

• We find strong infrared variability at all wave- 
lengths in every source. These fiuctuations are typ- 
ically rapid (days to weeks) and often show a wave- 
length dependence where the short-wavelength flux 
increases as the long- wavelength flux decreases. 

• The most likely source of this infrared variability is 
a vertical perturbation of the inner disk, based on 
the wavelength dependence and rapid timescales of 
the variability, as well as direct measurements of 
the emission from the inner disk. 



• The strength and timescale of the infrared variabil- 
ity in these evolved disks is not unlike that seen in 
the rest of the IC 348 cluster members with disks, of 
which roughly 60% are measurably variable. This 
result suggests that these types of perturbations 
are common. 

• It is unlikely that variable heating from fluctuations 
in the stellar flux or accretion rate are responsi- 
ble for the infrared variability. We do not observe 
strong stellar flux changes on the same timescale 
as the infrared variability, and the accretion lumi- 
nosity is not a significant component of the fiux 
striking the disk. 

• We can also exclude the influence of strong winds 
and companions at certain orbital radii within the 
disk. 

• Possible explanations for the behavior are either 
companions at permitted orbital radii, effects from 
the system magnetic field, or the influence of very 
large X-ray flares. However these models predic- 
tions that are not yet confirmed - in the first case 
of periodicities, in the second of the variable ex- 
tinction of the star and in the third of correlated 
X-ray and infrared flux. 

The size and timescale of the infrared flux variations 
seen in the six evolved disks studied here are very sim- 
ilar to the other variable disks in IC 348, which make 
up roughly half of the cluster members with disks. This 
suggests that many young stellar objects undergo similar 
perturbations to the structure of their inner disks. Fur- 
ther observations are needed to determine if the physical 
causes of the variability are similar between these two 
samples, but the similarity in the size of the fluctuations 
suggest that variability could be a useful tool for study- 
ing the inner disks around many young stellar objects. 
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TABLE 1 
Observing Log 



Date MJD Wavelength Resolution Instrument 



LRLL 2 


08 Oct 2007 


54381.9 


5-40//m 


~600 


IRS spectra 


15 Oct 2007 


54388.9 


5-40/im 


~600 


IRS spectra 


24 Mar 2008 


54549.2 


5-40/im 


~600 


IRS spectra 


31 Mar 2008 


54556.7 


5-40/im 


~600 


IRS spectra 


08 Oct 2009=^ 


55112.6 


0.8-2.5/im 


~800 


Spex 


31 Oct 2009^ 


55137.3 


0.8-2. 5/im 


~800 


Spex 


LRLL 21 


09 Oct 2007 


54382.3 


5-40ynm 


~600 


IRS spectra 


15 Oct 2007 


54388.8 


5-40/im 


~600 


IRS spectra 


24 Feb 2008 


54520.3 


5-40/im 


~600 


IRS spectra 


29 Feb 2008 


54525.1 


6460-6650A 


~ 35000 


Hectochelle 


01 Mar 2008 


54526.9 


5-40/im 


~600 


IR S spectra 


10 Oct 2008 


54750.5 


0.8-5//m 


~1500 


Spex 


11 Oct 2008 


54751.5 


0.8-2.5/im 


~ 800 


Spex 


18 Oct 2008 


54757.4 


O.S-bfim 


~1500 


Spex 


19 Oct 2008 


54758.4 


0.8-5/im 


~1500 


Spex 


08 Oct 2009=* 


55112.6 


0.8-5/im 


~ 800 


Spex 


12 Oct 2009=^ 


55116.1 


6000-9000A 


~ 3500 


CAFOS 


31 Oct 2009=^ 


55137.3 


0.8-5/im 


~ 800 


Spex 


04 Nov 2009^ 


55141.5 


0.8-5/im 


~ 800 


Spex 


08 Nov 2009 


55144.5 


0.8-2.5/im 


~ 800 


Spex 


LRLL 58 


08 Oct 2007 


54381.9 


5-40/im 


~600 


IRS spectra 


15 Oct 2007 


54388.9 


5-40/jm 


~600 


IRS spectra 


29 Feb 2008 


54525.1 


6460-6650A 


~ 35000 


Hectochelle 


24 Mar 2008 


54549.2 


5-40/im 


~600 


IRS spectra 


31 Mar 2008 


54556.8 


5-40/im 


~600 


IRS spectra 


09 Oct 2009=* 


55113.6 


0.8-2.5/tm 


~800 


Spex 


13 Oct 2009=^ 


55117.1 


6000-9000A 


~3500 


CAFOS 


08 Nov 2009 


55144.5 


0.8-2.5/im 


~800 


Spex 


LRLL 67 


09 Oct 2007 


54382.4 


5-40/im 


~600 


IRS spectra 


15 Oct 2007 


54388.9 


5-40/jm 


~600 


IRS spectra 


24 Feb 2008 


54520.3 


5-40/jm 


~600 


IRS spectra 


29 Feb 2008 


54525.1 


6460-6650A 


~ 35000 


Hectochelle 


03 Mar 2008 


54527.0 


5-40/tm 


~600 


IRS spectra 


LRLL 1679 


08 Oct 2007 


54381.9 


5-40/tm 


~600 


IRS spectra 


15 Oct 2007 


54388.9 


5-40/tm 


~600 


IRS spectra 


29 Feb 2008 


54525.1 


6460-6650A 


~ 35000 


Hectochelle 


24 Mar 2008 


54549.1 


5-40/im 


~600 


IRS spectra 


31 Mar 2008 


54556.8 


5-40/(m 


~600 


IRS spectra 


11 Oct 2008 


54751.5 


0.8-2.5/im 


~ 800 


Spex 


18 Oct 2008 


54757.4 


0.8-2.5/jm 


~1500 


Spex 


19 Oct 2008 


54758.4 


0.8-2.5/jm 


~1500 


Spex 


09 Oct 2009=* 


55113.6 


0.8-2.5/im 


~800 


Spex 


13 Oct 2009=^ 


55117.1 


6000-9000A 


~3500 


CAFOS 


04 Nov 2009=" 


55141.5 


0.8-2.5/tm 


~ 800 


Spex 


All Stars 


21 Feb 2004 


53056.6 


24/jm 


photometry 


MIPS 


19 Sep 2004 


53267.9 


24/im 


photometry 


MIPS 


23 Sep 2007 


54366.9 


24/im 


photometry 


MIPS 


24 Sep 2007 


54367.9 


24/im 


photometry 


MIPS 


25 Sep 2007 


54368.9 


24/jm 


photometry 


MIPS 


26 Sep 2007 


54369.9 


24/tm 


photometry 


MIPS 


27 Sep 2007 


54370.8 


24/(m 


photometry 


MIPS 


12 Mar 2008 


54538.7 


24/(m 


photometry 


MIPS 


19 Mar 2008 


54544.8 


24/tm 


photometry 


MIPS 


11 Feb 2004 


53046 


3-8/im 


photometry 


IRAC cold-mission 


08 Sep 2004 


53257 


3-8/tm 


photometry 


IRAC cold-mission 


19 Mar 2009 


54910 


3-8/im 


photometry 


IRAC cold-mission 


20 Mar 2009 


54911 


3-8/im 


photometry 


IRAC cold-mission 


21 Mar 2009 


54912 


3-8/im 


photometry 


IRAC cold-mission 


22 Mar 2009 


54913 


3-8/im 


photometry 


IRAC cold-mission 


23 Mar 2009 


54914 


3-8/im 


photometry 


IRAC cold-mission 


03 Oct-07 Nov 2009 


55107-55142 


3.6,4.5/tm 


photometry 


IRAC warm-mission 



^ Data were taken during the IRAC warm-mission monitoring campaign in fall 2009. 
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TABLE 2 

MIPS PHOTOMETRY 



MJD LRLL 2 LRLL 21 LRLL 58 LRLL 67 LRLL 1679 



53055.56 


3.288 


± 


0, 


,015 


3, 


,651 


± 


0, 


,021 


4.686 


± 


0.015 


4, 


,658 


± 


0.011 


6, 


,561 


± 


0, 


,009 


53266.87 


3.611 


± 


0, 


,019 


3, 


,866 


d= 


0, 


,014 


4.602 


± 


0.015 


4, 


,654 


± 


0.014 


7, 


,025 


± 


0, 


,013 


53267.09 


3.587 




0, 


,017 


3, 


,872 


± 


0, 


,013 


4.557 


± 


0.010 


4, 


,644 


± 


0.017 


6, 


,945 


± 


0, 


,014 


54365.88 


3.413 


± 


0, 


,013 


3, 


,717 


± 


0, 


,011 


4.531 


± 


0.011 


4, 


,649 


± 


0.012 


6, 


,972 




0, 


,012 


54366.86 


3.407 


± 


0, 


,020 


3, 


,753 


d= 


0, 


,011 


4.528 


± 


0.014 


4, 


,688 




0.012 


6, 


,925 


± 


0, 


,014 


54367.91 


3.409 


± 


0, 


,019 


3, 


,741 


± 


0, 


,012 


4.520 


± 


0.012 


4, 


,668 




0.013 


6, 


,893 


± 


0, 


,012 


54368.86 


3.402 


± 


0, 


,013 


3, 


,719 


± 


0, 


,010 


4.526 


± 


0.013 


4, 


,667 


± 


0.011 


6, 


,953 


± 


0, 


,011 


54369.80 


3.386 


± 


0, 


,016 


3, 


,778 


± 


0, 


,007 


4.546 


± 


0.017 


4, 


,649 


± 


0.016 


7, 


,128 


± 


0, 


,013 


54537.66 


3.367 


± 


0. 


,013 


3, 


,631 


± 


0, 


,020 


4.581 


± 


0.016 


4, 


.656 


± 


0.010 


6, 


.563 


± 


0, 


Oil 


54543.78 


3.337 


± 


0. 


,015 


3, 


,730 


± 


0, 


,016 


4.546 


± 


0.014 


4, 


.650 


± 


0.012 


6, 


.631 


± 


0, 


Oil 



TABLE 3 
IHAC I'lioi OMKTin" 



Date [3.6] [4.5] [5.8] [8.0] 



11 Feb 2004 
08 Sep 2004 

19 Mar 2009 

20 Mar 2009 

21 Mar 2009 

22 Mar 2009 

23 Mar 2009 



6.971 ± 0.002 

7.015 ± 0.002 
6.982 ± 0.002 
6.995 ± 0.002 
6.948 ± 0.002 
6.984 ± 0.002 



LRLL 2 



6.797 ± 0.002 

6.905 ± 0.003 
6.791 ± 0.002 
6.748 ± 0.002 
6.796 ± 0.002 
6.800 ± 0.002 



6.555 ± 0.002 

6.455 ± 0.002 
6.454 ± 0.002 
6.479 ± 0.002 
6.468 ± 0.002 
6.500 ± 0.002 



5.761 ± 0.004 

5.724 ± 0.005 
5.721 ± 0.005 
5.687 ± 0.004 
5.713 ± 0.003 
5.742 ± 0.004 



LRLL 21 



11 Feb 2004 


9.086 


± 


0.003 


8.912 


it 


0.003 


8.619 


it 


0.006 


8.130 


it 


0.02 


08 Sep 2004 


9.069 


± 


0.004 


8.657 


± 


0.003 


8.413 


± 


0.006 


8.014 


it 


0.02 


19 Mar 2009 


9.221 


± 


0.003 


9.025 


± 


0.003 


8.862 


± 


0.007 


8.306 


it 


0.03 


20 Mar 2009 


9.212 


± 


0.003 


9.033 


± 


0.003 


8.806 


± 


0.006 


8.313 


it 


0.03 


21 Mar 2009 


9.119 


± 


0.003 


8.961 


± 


0.003 


8.787 


± 


0.006 


8.270 


it 


0.03 


22 Mar 2009 


9.063 


± 


0.003 


8.885 


± 


0.003 


8.669 


± 


0.005 


8.193 


it 


0.02 


23 Mar 2009 


8.934 


± 


0.003 


8.746 


± 


0.003 


8.511 


± 


0.005 


8.168 


it 


0.02 



LRLL 58 



11 Feb 2004 


9.789 


it 


0.004 


9.588 


± 


0.004 


9.330 


it 


0.015 


8.543 


it 


0.05 


08 Sep 2004 


10.104 


: it 


0.004 


9.959 


it 


0.004 


9.637 


it 


0.017 


8.557 


it 


0.05 


19 Mar 2009 


9.605 


it 


0.004 


9.372 


it 


0.003 


9.059 


it 


0.014 


8.513 


it 


0.04 


20 Mar 2009 


9.619 


it 


0.003 


9.347 


± 


0.003 


9.043 


it 


0.012 


8.493 


it 


0.04 


21 Mar 2009 


9.703 


it 


0.004 


9.305 


± 


0.003 


9.011 


± 


0.014 


8.506 


it 


0.05 


22 Mar 2009 


9.556 


it 


0.004 


9.254 


± 


0.004 


8.973 


± 


0.012 


8.484 


it 


0.05 


23 Mar 2009 


9.515 


it 


0.003 


9.245 


± 


0.004 


8.939 


it 


0.012 


8.494 


it 


0.05 



LRLL 67 



11 Feb 2004 


10.425 


it 


0.003 


10.255 


it 


0.003 


10.063 


it 


0.007 


9.575 it 0.010 


08 Sep 2004 


10.425 


it 


0.003 


10.238 


it 


0.004 


9.986 


it 


0.007 


9.509 it 0.011 


19 Mar 2009 


10.422 


it 


0.003 


10.343 


it 


0.003 


10.225 


it 


0.008 


9.675 it 0.010 


20 Mar 2009 


10.325 


it 


0.002 


10.284 


it 


0.003 


10.129 


it 


0.008 


9.633 it 0.010 


21 Mar 2009 


10.524 


it 


0.003 


10.325 


it 


0.003 


10.149 


it 


0.008 


9.644 it 0.010 


22 Mar 2009 


10.386 


it 


0.003 


10.439 


it 


0.004 


10.111 


it 


0.008 


9.605 it 0.009 


23 Mar 2009 


10.426 


it 


0.003 


10.388 


it 


0.004 


10.163 


it 


0.008 


9.658 it 0.012 










LRLL 


1679 










11 Feb 2004 


11.317 


it 


0.004 


11.060 


it 


0.004 


10.781 


it 


0.011 


10.287 it 0.010 


08 Sep 2004 


10.806 


it 


0.003 


10.446 


it 


0.003 


10.173 


it 


0.008 


9.809 it 0.007 


19 Mar 2009 


11.468 


it 


0.004 


11.295 


± 


0.005 


11.063 


± 


0.012 


10.769 it 0.014 


20 Mar 2009 


11.436 


± 


0.004 


11.267 


± 


0.005 


11.099 


± 


0.013 


10.781 it 0.014 


21 Mar 2009 


11.459 


± 


0.004 


11.305 


it 


0.005 


11.047 


it 


0.012 


10.744 it 0.015 


22 Mar 2009 


11.449 


it 


0.004 


11.298 


± 


0.005 


11.074 


± 


0.013 


10.728 it 0.016 


23 Mar 2009 


11.444 


it 


0.004 


11.303 


it 


0.005 


11.093 


it 


0.013 


10.769 it 0.015 
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TABLE 4 
IRAC Warm Mission Photometry 



LRLL 2 


LRLL 21 


LRLL 58 


LRLL 67 


LRLL 1679 


MJD-55000 [3.6] [4.5] 


[3.6] [4.5] 


[3.6] [4.5] 


[3.6] [4.5] 


[3.6] [4.5] 



lUD. 


Q1 

.y 1 


7 


01 

.UIO 


o. 


Sfi9 
.oOZ 


S OQfi 

o.yyo 


S SOQ 

o.ouy 




y. 


7A^ 




y. 


4SQ 

.^oy 


1 

lU. 


4S7 
.^o i 


1 

lU. 


401 
,^U1 


1 1 

1 1 , 


97S 
.Z / o 


1 1 
1 1 


09Q 

.uzy 


1 07 


. ou 


7 

1 . 


OOQ 

. uuy 


fi 

U. 


. ouo 


8 979 


8 804 


Q 

y . 


77S 
■ 1 i o 


Q 

y . 


4Qfi 

.^yu 


1 

lU. 


•^00 

.ouu 


1 
lU. 


41 7 


11 


.OZO 


11 


090 

.uzu 


1 07 

J-U 1 ■ 
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7 

1 . 
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fi 

u. 


. ooo 
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y .uz i 
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y . 


4fiQ 


1 

lU. 


.ouo 
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lU. 
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11 
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1 
lU 


QOQ 

.yyy 


1 07 
lU ( . 


.ol 


fi 

D. 
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fi 

u. 


S70 


Q OOfi 

y .uuD 
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y. 
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1 

lU. 


f^91 
.OZl 


1 

lU. 
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. lo 
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1 . 
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o.yy 1 
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Q 

y. 
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Q 

y. 


.400 


1 
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477 
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S Q71 
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Q 
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794 
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Q 
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1 
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1 
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,4:00 




.04:0 




0Q9 

.uyz 


1 OS 


.Oo 


O. 


QQ1 

.yy 1 
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Q Ci^A 
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Q 
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Q 
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Q 
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1 
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1 
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Q4S 
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TABLE 5 
Extinction 



Date 






LRLL 2 


08 Oct 2009 


2.8 


31 Oct 2009 


3.0 


LRLL 21 


10 Oct 2008 


4.1 


11 Oct 2008 


4.4 


18 Oct 2008 


4.2 


19 Oct 2008 


4.6 


08 Oct 2009 


3.9 


31 Oct 2009 


4.1 


04 Nov 2009 


4.2 


08 Nov 2009 


4.3 


LRLL 58 


09 Oct 2009 


3.4 


08 Nov 2009 


3.4 


LRLL 1679 


11 Oct 2008 


5.8 


18 Oct 2008 


5.8 


19 Oct 2008 


5.8 


09 Oct 2009 


5.8 


04 Nov 2009 


5.8 



TABLE 6 

IlMISSIOX LiXES 



Date Pa/3 EW (A) Br7 EW (A) Corrected Br7 EW (A) 



LRLL 2 


9 Oct 2009 
31 Oct 2009 


-0.46 ± 0.04 
-3.40 ± 0.03 


-3.45 ± 0.27 
-5.23 ± 0.17 


-1.85 ± 0.17 
-3.16 ± 0.10 


LRLL 21 


10 Oct 2008 

11 Oct 2008 

18 Oct 2008 

19 Oct 2008 
8 Oct 2009 
31 Oct 2009 
4 Nov 2009 
8 Nov 2009 


0.63 ± 0.35 
0.96 ± 0.07 
-1.21 ± 0.05 
-1.25 ± 0.08 
0.83 ± 0.06 
1.90 ± 0.07 
-0.70 ± 0.07 
-0.56 ± 0.06 


-0.01 ± 0.27 
0.32 ± 0.05 
0.16 ± 0.04 
0.07 ± 0.05 
-0.29 ± 0.05 
-0.14 ± 0.04 
-1.02 ± 0.08 
-0.88 ±0.09 


-0.07 ± 0.26 
0.27 ± 0.04 
0.09 ± 0.03 
0.02 ± 0.04 
-0.33 ± 0.05 
-0.22 ± 0.03 
-1.11 ± 0.08 
-0.94 ± 0.09 


LRLL 58 


9 Oct 2009 
8 Nov 2009 


-0.96 ± 0.09 
-1.95 ± 0.51 


-0.66 ± 0.09 
-1.63 ± 0.45 


-0.86 ± 0.09 
-1.84 ± 0.47 


LRLL 1679 


11 Oct 2008 

18 Oct 2008 

19 Oct 2008 
9 Oct 2009 
4 Nov 2009 


0.58 ± 0.22 
0.23 ± 0.24 
1.03 ± 0.33 
0.89 ± 0.49 
0.78 ± 0.31 


0.12 ± 0.07 
-0.66 ± 0.20 
-0.13 ± 0.15 
-0.51 ± 0.42 
-0.17 ± 0.13 





TABLE 7 
Accretion Rates 



Date 


Flux (10-^"' erg s"^ cm-^) 


log(L/L0) 


log(Lacc/L0) 




M (10-«MQj/r-i) 




LRLL 2 Pa/3 


9 Oct 2009 


25.2 ± 6.8 


-3.09 ± 0.12 


-0.37 ± 0.14 


0.003 


3.14 ± 0.44 


31 Oct 2009 


185 ± 29 


-2.23 ± 0.07 


0.61 ± 0.08 


0.03 


29.9 ± 2.4 


LRLL 2 Br7 


9 Oct 2009 


29.2 ± 3.6 


-3.03 ± 0.05 


0.61 ± 0.06 


0.03 


29.9 ± 1.8 


31 Oct 2009 


49.8 ± 5.2 


-2.80 ± 0.04 


0.90 ± 0.05 


0.06 


58.4 ± 2.9 


LRLL 58 Pa/3 


9 Oct 2009 


1.35 ± 0.24 


-4.37 ± 0.08 


-1.83 ± 0.09 


0.02 


0.23 ± 0.02 


8 Nov 2009 


2.81 ± 0.56 


-4.05 ± 0.09 


-1.47 ± 0.10 


0.05 


0.52 ± 0.05 


LHLL o,s Bi- 


9 Oct 2009 


0.74 ± 0.09 


-4.63 ± 0.06 


-1.40 ± 0.08 


0.06 


0.62 ± 0.05 


8 Nov 2009 


1.57 ± 0.23 


-4.30 ± 0.06 


-0.99 ± 0.08 


0.14 


1.58 ± 0.13 


LRLL 21 Pa/3 


10 Oct 2008 


< 3.52 


< -3.95 


< -1.35 


< 0.02 


< 0.26 


11 Oct 2008 


< 0.77 


< -4.61 


< -2.11 


< 0.003 


< 0.05 


18 Oct 2008 


4.68 ± 0.74 


-3.83 ± 0.07 


-1.22 ± 0.08 


0.02 


0.388 ± 0.031 


19 Oct 2009 


5.33 ± 0.86 


-3.77 ± 0.07 


-1.15 ± 0.08 


0.02 


0.455 ± 0.036 


8 Oct 2009 


< 0.79 


< -4.60 


< -2.09 


< 0.002 


< 0.05 


31 Oct 2009 


< 0.95 


< -4.52 


< -2.01 


< 0.003 


< 0.07 


4 Nov 2009 


3.36 ± 0.57 


-3.97 ± 0.07 


-1.38 ± 0.09 


0.01 


0.268 ± 0.024 


8 Nov 2009 


2.75 ± 0.49 


-4.06 ± 0.08 


-1.48 ± 0.10 


0.01 


0.212 ± 0.021 


LRLL 21 Br7 


10 Oct 2008 


< 1.14 


< -4.34 


< -1.04 


< 0.003 


< 0.53 


11 Oct 2008 


< 0.22 


< -5.15 


< -2.06 


< 0.002 


< 0.05 


18 Oct 2008 


< 0.19 


< -5.23 


< -2.16 


< 0.003 


< 0.05 


19 Oct 2008 


< 0.26 


< -5.09 


< -1.98 


< 0.04 


< 0.07 


8 Oct 2009 


0.663 ± 0.164 


-4.68 ± 0.11 


-1.47 ± 0.14 


0.01 


0.218 ± 0.031 


31 Oct 2009 


0.501 ± 0.103 


-4.79 ± 0.09 


-1.61 ± 0.11 


0.007 


0.158 ± 0.017 


4 Nov 2009 


2.78 ± 0.31 


-4.05 ± 0.05 


-0.67 ± 0.06 


0.06 


1.38 ± 0.08 


8 Nov 2008 


2.30 ± 0.29 


-4.13 ± 0.05 


-0.77 ± 0.06 


0.05 


1.09 ± 0.07 


LRLL 1679 Pa/3 


11 Oct 2008 


< 0.35 


< -4.96 


< -2.49 


< 0.02 


< 0.08 


18 Oct 2008 


< 0.38 


< -4.92 


< -2.46 


< 0.02 


< 0.08 


19 Oct 2008 


< 0.52 


< -4.78 


< -2.30 


< 0.02 


< 0.12 


9 Oct 2009 


< 0.78 


< -4.61 


< -2.10 


< 0.04 


< 0.18 


4 Nov 2009 


< 0.49 


< -4.81 


< -2.33 


< 0.07 


< 0.11 


LRLL 1679 Br7 


11 Oct 2008 


< 0.06 


< -5.70 


< -2.75 


< 0.008 


< 0.04 


18 Oct 2008 


< 0.18 


< -5.25 


< -2.19 


< 0.03 


< 0.15 


19 Oct 2008 


< 0.13 


< -5.37 


< -2.34 


< 0.02 


< 0.11 


9 Oct 2009 


< 0.38 


< -5.92 


< -3.03 


< 0.004 


< 0.02 


4 Nov 2009 


< 0.12 


< -5.43 


< -2.41 


< 0.02 


< 0.09 



TABLE 8 
LRLL 21 NIR Photometry 



Date 


J 


H 


K 


2MASS 


11.02 


9.99 


9.47 


10 Oct 2008 


11.31 


10.17 


9.71 


11 Oct 2008 


11.30 


10.21 


9.73 


19 Oct 2008 


11.19 


10.14 


9.60 


8 Oct 2009 


10.96 


9.98 


9.59 


31 Oct 2009 


10.99 


10.05 


9.52 


4 Nov 2009 


10.95 


9.91 


9.39 


8 Nov 2009 


10.96 


9.94 


9.43 
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TABLE 9 
LRLL 21 Stellar 
Luminosity 



Date L* (L©) 



10 Oct 2008 

11 Oct 2008 
19 Oct 2008 
08 Oct 2009 
31 Oct 2009 
04 Nov 2009 
08 Nov 2009 



2.52 pm 0.39 
2.75 pm 0.43 
3.20 pm 0.50 
3.30 pm 0.52 
3.38 pm 0.53 
3.60 pm 0.56 
3.66 pm 0.57 



TABLE 10 
Radial and Rotational Velocity 



ID Vr vsini Reference 



LRLL 21 


17.01 


17.84 


1 




-15.7=^ 


21.9'' 


2 




66. 8=^ 


<15'' 


2 


LRLL 58 


14.2 


<15 


2 


LRLL 67 


15.0 


<15 


2 




15.09 


13.6 


1 




16.9 


<11 


3 



Not e. — Referen ces: (1) iDah^ 
I 200^ ) (2) This paper O^ INordhagen et al.| 
\ 20m) . 

^ In our high-resolution spectra of LRLL 
21 we find two peaks in the cross- 
correlation function and report the radial 
and rotational velocity of each component. 
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TABLE 11 
Infrared Excess 



Date 










Covering Fraction 




LRLL 21 


10 Oct 2008 


0.13 


0.19 


0.26 


1900 


0.06 


11 Oct 2008 


0.10 


0.22 


0.35 






18 Oct 2008 


0.01 


0.17 


0.36 


2020 


0.08 


19 Oct 2008 


0.01 


0.17 


0.37 


2180 


0.08 


08 Oct 2009 


0.02 


0.09 


0.19 




~ 0.04^^ 


31 Oct 2009 


0.06 


0.29 


0.56 


1860 


0.13 


04 Nov 2009 


0.07 


0.29 


0.54 


1700 


0.13 


08 Nov 2009 


0.15 


0.29 


0.46 






LRLL 58'' 


09 Oct 2009 


0.12 


0.28 


0.46 




~ 0.13=^ 


08 Nov 2009 


0.18 


0.32 


0.48 







LRLL 1679 



11 Oct 2008 


-0.02 


-0.004 


0.01 






18 Oct 2008 


-0.09 


-0.07 


-0.04 






19 Oct 2008 


-0.16 


-0.11 


-0.05 






09 Oct 2009 


-0.11 


-0.02 


0.07 




~ 0.02=^ 


04 Nov 2009 


-0.04 


0.02 


0.09 




~ 0.03a 



Note. — Uncertainties on the veiling measurements are ~0.1. No 
veiling is measured for LRLL 2 since the only photospheric absorption 
lines are the hydrogen series, which will be contaminated from emission 
by the accretion flow. 

^ For days when 2-5jU.m spectra were not available, but we have Spitzer 
3.6 and 4.5/.im photometry, we estimate a blackbody fit to the excess 
assuming T=1900K, a temperature that was found to fit those days 
that do have the 2-5jum spectra. 

^ For LRLL 58 we use the standard JH 108 to measure the veiling 
even though it shows emission near the Br7 line, which may indicate 
accretion and possibly a disk. The line strengths fit better than an Ml 
giant or dwarf standard. The measured veiling is only a lower limit, 
although the lack of infrared excess for JH 108 indicates that within 
our uncertainties we are measuring the true veiling 
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TABLE 12 
Variability Summary 



Star 


Stellar Flux 




Accretion Flux 




Infrared Variability 


(Spectral Type) 


Variable?'' 


# of Epochs 


Variable? Rate 


(lO-*M0yr-l) 


Strength 


Pivot Point 




LRLL 2 
A2 


no 


2 


no?'^ 


~30 


20% in IRS 
25% at [24] 
0.1 mag at [3. 6], [4.5] 


6/^m 


LRLL 21 
KO 


P: 2.5 days 
NP: yearly 


8 


yes 


< 0.05 - 0.46 


40% for A < 14/xm in IRS 
< 8% for A > 14Atm in IRS 
zU/o at [z4J 
0.5 mag at [3. 6], [4.5] 


8/im 


LRLL 31 
G6 


P: 3.4 days 


9 


yes 


0.25-1.62 


60% in IRS 
30% at [24] 
0.3 mag at [3. 6], [4. 5] 


8/im 


LRLL 58 
Ml. 25 


P: 7.4 days 


2 


yes 


0.23-0.52 


<5% for A < 14/xm in IRS 
< 20% for A > Ufj.m in IRS 
14% at [24] 
0.3 mag at [3. 6], [4.5] 


8/im 


LRLL 67 
MO. 75 


no 


3^ 


yes^ 


~0.2<^ 


< 3% for A < 14/.tm in IRS 

< 20% for A > Mfim in IRS 

40% at [24] <= 
0.4 mag at [3.6], [4.5] 


none detected 


LRLL 1679 
M3.5 


unknown 


5 


undetected 


< 0.18 


50% in IRS 
50% at [24] 
0.1 mag at [3.6], [4.5] 


13/xm 



^ P— Periodic fluctuations. NP— fluctuations detected with no known period. For variable sources we include the period or the typical timescale of 
the non-periodic variability. 

^ Derived from our observations of the Pa/3 line. 

The large systematic uncertainties in characterizing the continuum around the Pa/3 line prevent us from determining if the accretion rate is variable 
even though the derived accretion rates in the two epochs are substantially different. 

^ We have no Pa/3 spectra for LRLL 67, which prevents us from measuring an accretion rate. The accretion rate is taken from [P^ahml <]2008l ) and the 

presence of variability is based on the changing Ha EW. 

^ The 24/j.m variability of LRLL 67 is only detected at the ~2cr level 
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TABLE 13 

Strength and shape of silicate feature 



Epoch EW(lOAjm) S(peak) S{11.3)/5(9.8) 



LRLL 2 



1 


2.98 


0.97 


1.00 


2 


2.71 


0.91 


1.06 


3 


2.97 


0.96 


1.11 


4 


2.96 


0.96 


1.13 


LRLL 21 


1 


2.19 


0.76 


1.39 


2 


2.13 


0.80 


1.52 


3 


1.87 


0.65 


1.53 


4 


1.89 


0.68 


1.81 


LRLL 31 


1 


2.37 


1.98 


0.98 


2 


3.31 


2.26 


0.96 


3 


2.37 


1.84 


0.97 


4 


2.37 


1.71 


0.99 


LRLL 58 


1 


2.51 


0.79 


1.13 


2 


2.76 


0.87 


1.18 


3 


3.07 


1.03 


1.27 


4 


3.21 


1.05 


1.24 


LRLL 67 


1 


3.73 


1.62 


0.87 


2 


3.51 


1.49 


0.87 


3 


3.92 


1.72 


0.87 


4 


3.77 


1.63 


0.86 


LRLL 1679 


1 


1.68 


0.69 


0.93 


2 


1.33 


0.47 


1.16 


3 


1.74 


0.69 


0.70 


4 


1.19 


0.53 


1.21 



Note. — These parameters have been derived based 
on a dcrcddened silicate feature. The extinction is as- 
sumed to be constant betwe en each epoch an d the spec- 
tra arc dcrcddened using the lMcCluTg 1 120091) extinction 
law 



TABLE 14 
Inner radii of disk 



ID Co-rotation Radius (AU) Sublimation Radius (AU) 

LRLL 2 (0.03-0.18) 0.76 (>0.41) 

LRLL 21 0.04 0.13 {>0.07) 

LRLL 31 0.05 0.13 (>0.07) 

LRLL 58 0.07 0.05 {>0.03) 

LRLL 67 (0.02-0.10) 0.04 (>0.02) 

LRLL 1679 (0.01-0.08) 0.03 (>0.02) 



Note. — Co-rotation radii arc calculated based on observed rotational 
periods estimated from optical light curves. When no observed period is 
available we estimate the co-rotation radius for periods of 1 and 15 days, 
which span tlic typical periods observed for prc-main sequence stars in IC 
348 (ICieza & Baliberii2006iV The lower limit on the sublimation radius, 
included in parenthesis, assumes very large grains (> lO/^m) at 1500K. 



